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The decay constants (k ;) of the equation of inert gas exchanges are 


the roots of an algebraic equation of degree n + 1, where n is the num- 
ber of distinct absorbing tissues. The coefficients of this equation can 
be obtained numerically by certain independent experiments to measure 
the tissue parameters. Graphical solution of this equation yields theo- 
retical values of the k;. Combining these constants with the numerical 


values for the partial derivatives of the k, then gives the per cent rate 
of change of the k; as any one tissue parameter varies by a given frac- 


tion of its normal range. A numerical example of these calculations 
shows good conformity with experiment, and permits a quantitative esti- 
mate of variations in the speed of gas exchange from a knowledge of 
changes in the physiological state. 


A central problem in aereo- and submarine physiology is the con- 
trol of the rate of inert gas exchange in such a way that personnel 
do not develop the “bends” or “chokes”. For either selection or thera- 
peutic purposes it is therefore indispensable to know the factors 
which govern the exchange and the effectiveness of each. The fol- 
lowing is an attempted solution of the problem based on an elemen- | 
tary mathematical analysis which we have presented elsewhere 
(Smith and Morales, 1944a, 1944b; Morales and Smith, 1944). 

Tissue regions (Figure 1), of which the entire body is a distinct 
parallel arrangement,? absorb? inert gases in a manner described by 


the equation: 


* The opinions expressed in this article are the private ones of the writers, 
and are not to be construed as reflecting the policies of the Navy Department or 
the Naval Service at large. b 

1 We have shown elsewhere (Morales and Smith, 1945) that tissues may be 
arranged in three essentially distinct ways, of which this is one. The other two 
arrangements predict the loss of two exponentials at a time when the number of 
distinct absorbing tissues is reduced by one. Under similar circumstances, with 
a competitive parallel arrangement, only one exponential should be lost. As A. 
Behnke pointed out in personal communication, the use of helium is practically 
tantamount to reducing the number of distinct tissues by one, and helium ex- 
change can be described by one less exponential than the exchange of nitrogen or 
krypton. This fact, plus the evidence alluded to later, forms the basis for our 


choice of arrangement. *: Ey 
2 Or, with suitable changes in the boundary conditions, desorbs. 
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¢(t) = (0) — OQ; 6 — G5 ie et 7 Minn 9 PY (1) 


where ¢(t) is the amount of gas taken up to time t, and ¢(c) is 
the asymptotic (steady state) amount; ” is the number of tissues, 
and the Q’s and the k’s are constants depending on the physiological 
state of the limb and the properties of the gas. 


SCHEMA OF A TISSUE REGION 
COMPETITIVE PARALLEL 


PERMEABILITY OF TOTAL BARRIER, 
BLOOO-FAT «ho 


TOTAL CAPILLARY SURFACE OF FAT = $5 
TOTAL VOLUME OF FAT= V2 


OF TOTAL BARRER, 
BLOOD-MUSCLE = hy 


TOTAL CAPILLARY SURFACE OF MUSCLE: $s, (2) 
TOTAL VOLUME OF MUSCLE = v, %, 
TRUE VOLUME OF MUSCLE = V, 


RATE OF BLOOD FLOW: R OF BLOOD FLOW? R 
ARTERIAL CONCENTRATION © C VENOUS CONCENTRATION © Xq 


FIGURE 1 


We consider that the parameters of the phenomenon are the fol- 
lowing quantities: 


Rem’ sec] = rate of blood flow through the limb, 


Vo[em*] = total amount of fluid outside of cells, including the 

blood volume, 
and for the ith tissue: 
a; = concentration in blood/concentration within cells, 
at equilibrium. 
h,[em sec] = cell membrane permeability, 
V;i[em*] = total cellular volume, 
Si 


v. [cm-?] = capillary surface/gross volume of tissue. 
ry 


It is furthermore convenient to define the function, 


Si 
nay — aj h; ee 

Now it has been shown by several authors (Underwood and Diaz, 
1941; Jones, 1941; Smith and Morales, 1944b) that equation (1) tak- 
en as a purely empirical fit, describes gas exchange data with remark- 
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able precision. This is true in particular of measurements of radio- 
active inert gas exchange by various tissue regions (Lawrence-Jones 
technique), to which we will turn presently. In such a procedure the 
decay constants, k; , can be most conveniently determined by graphical 
analysis on arith-log paper. The order of magnitude of such constants 
is therefore well established. On the other hand, the parameters list- 
ed above can all be measured or estimated independently of gas ab- 
sorption experiments (vide infra). When the values of these para- 
meters are incorporated into the equations, the theoretical values of 
the k; may be computed. It is proposed here to demonstrate the theo- 
retical determination of decay constants, to show a favorable com- 
parison with empirically obtained values, and finally to consider the 
effect of varying the physiological state on the decay constants, hence 
on the speed of inert gas exchange. 

Values of the &; are determined by [Smith and Morales, 1944a, 
equation (9)]: 


heh h; S; k; 
Vous > 


is1 Ti — Wj 


—R=0, (2) 


the left member of which is a pelynomial of degree n + 1. In mam- 
malian systems the number (n) of absorbing tissues is practically 
two, i.e., a tissue region is composed of aqueous tissue, fat, and bone, 
the last being essentially non-absorbing. 

Clearing equation (2) of fractions for the special case, n = 2, 
we obtain 


Vi V2 
(+ a) Ve ta te TR 
Pe paetsected Serene Fo 
Vo 
vv (3) 
mam( Vor +?) +R(m +m) 
“5 Oy Oe (oe = 
v, j v. , 


where the subscripts 1 and 2 denote aqueous tissue and fat respec- 
tively. On substitution of values for the tissue parameters equation 
(3) becomes a numerical cubic, which on graphical solution gives 
three decay constants, k,, ki, ke. 

We may now differentiate equation (3) partially and obtain, 


102 INERT GAS EXCHANGE 


ker.) aie 
lok; = wo; k; 
eyo (Ry Sy). oUt leg ae 
10k; sf Os -\? 
= wy Ve we lie: (4) 
kj 0 as tg eye 
1d k; 
Ta OFs 
k; 0 Vo 
lok; 2 ( Ts y 
k;oVs Og tps; : 
with 
ia; a 
o; = im V, = : A) 
con ) 
ix1 Qi mi — ly 


corresponding to the per cent changes in the k; for cgs unit change 
in a tissue parameter. These derivatives need only be corrected in a 
physiological sense in order to answer fully the problem we proposed 
to solve. By this we mean, for example, that a change of 1 cm’ in 
volume of muscle is manifestly not the physiological equivalent of a 
change of 1 cm® sec! in rate of blood flow. We may, however, make 
these changes comparable by finding the per cent change in the k;,’s 
when tissue paramters vary by, say, 10% of their normal range. Mul- 
tiplying the differential coefficients by such factors then gives the 
final desired result. 


The methods for obtaining values of the tissue parameters are 
various, and will be considered in order. 

The average rate of blood flow over a three hour period can of 
course be measured by a water or air plethysmograph.° 

The total amount of fluid which is outside of the cells is the sum 
of three contributions: (1) The true blood volwme, which can be ap- 
proximated as being equal to the volume change of a limb in pass- 
ing from a very cold to a warm environment. (2) The free-space of 
muscle. Sinee Cl does not penetrate muscle cells, any “chloride 
space” method (Fenn, 1936) gives the free space fraction for muscle. 
This value (18%) can also be checked theoretically, assuming quasi- 
close packing of cylindrical muscle fibers. The free space fraction 

3 The egs dimensions for flow are cms® sec-!, but in mammalian physiology 


it is customary to express blood flow in cubie centimeters of blood per 100 cubic 
centimeters of tissue. 
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times the total volume then gives the free space. (38) The free space 
of fat. Since the chloride method has not been shown to be applicable 
to fat tissue for the reason that, unlike the situation in muscle, the 
Cl" may be absorbed by the cells, recourse was sought in the value 
for the theoretical free space fraction of close packed ellipsoids, which 
can be shown to be 26%.‘ Again, this fraction times the total volume 
yields total free space. 


The partition coefficients can be measured directly (Loomis, 
1941). 


The ratio of capillary surface/volume of tissue was judged from 
the work of I. Gersh and M. A. Still (1945), with due allowance for 
the differences which can be expected to exist in the human hand. 

The plasma membrane permeabilities probably involve the great- 
est error of all of these determinations, since for virtually all systems 
they are known little more than by the order of magnitude. The com- 
putations of H. D. Landahl (1939) show oxygen permeability to be 
of the order of 10° cm sec. Furthermore, the data of F. M. Miiller 
(1941) and of others (see International Critical Tables) show that 
the permeabilities of oxygen and krypton are probably very similar. 
Guided by Professor S. C. Brooks’ opinion that fat cells may have a 
higher h for krypton than do muscle cells, we have tentatively as- 
sumed 2 X 10° and 1 X 10° respectively as the krypton plasma mem- 
brane permeabilities of fat tissue and muscle. 

It will be noted that when dealing with surface, we have spoken 
about capillary surface, while in dealing with permeability we have 
used plasma membrane values. This we feel is justified because the 
physiological behavior of the gross barrier will be governed by the 
limiting factors which go to make it up.’ If such is the case, then 
certainly the limiting surface is that of the capillaries, and the limit- 
ing permeability that of the plasma membrane. 

The total volumes were estimated from roentgenographs of 
hands. 

Use of the foregoing techniques can be illustrated by a case of 
absorption of radioactive krypton by the hand. When measurements 


4 This is based on a simple generalization of the calculation which leads to 
the same values for close-packed spheres. In histological section, fat cells appear 
compressed, hence not subject to being treated as rigid ellipsoids; however, a 
number of observations by Ens. Mary Still and one of us (MFM) were enough 
to show that surviving fat cells in isotonic solution are truly arranged in the 
hexagonal lattice of rigid ellipsoids. 

5 This can be easily shown analytically by solving the transfer problem for 
a double barrier. If the permeability of the first membrane and the surface of 
the second are assumed to be much greater than the corresponding quantities for 
the other membrane, then the problem degenerates to one of transfer across a 
membrane having the lower permeability and the smaller surface. 
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VALUES ON WHICH CALCULATIONS ARE BASED: 


STANDARD C.G.S. EQUIVALENT OF 
PARAMETER VALUE RANGE 10% PHYSIOLOGICAL RANGE 
R RATE OF BLOOD FLOW 13 10.0 1.0 
Vo BLOOD CHAMBER VOLUME 100.0 100.0 10.0 
hs, MUSCLE PENETRATION 0.12 1.188 0.1188 
100 Y, MUSCLE PARTITION 1.0 0.2 0.02 
Vv, MUSCLE VOLUME 108.2 100.0 10.0 
hySp FAT PENETRATION 0.044 0.4356 0.04356 
90 %. FAT PARTITION 0.125 0.95 0.095 


FAT VOLUME 10.0 


STANDARD K 
Kot.I7_ min"! 
K, =.0285 

Ke #0029 


Ses 50 


: k 
of *2 
-10 
(Vp) (m4) (4) (y) mt? (Gy) (Vp) 
RATE BLOOD MUSCLE MUSCLE MUSCLE FAT FAT 
-20 Eos CHAMBER PENETRATION PARTITION VOLUME PENETAATION PARTITION =: VOLUME 
FLOW [cms3]  [ems3sec-') [ers3] — [crms3sec?] [ems] 
FIGURE 2 


of subject C.J.S. (tabulated in Figure 2) are substituted into equa- 
tion (3), we obtain, 


kj? — (3.87582 X 10+) k,? 
+ (1.519854 X 10-°) k; — (.06643 X 10°) =0, 


whose graphical solution appears in Figure 3. On the other hand, 
when we substitute the values into the partial derivatives (4) and 


GRAPHICAL SOLUTION OF THE AUXILIARY EQUATION 


+1 


a 


FIGURE 3 
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correct for range, we obtain the per cent changes (Figure 2) induced 
in the k; (hence on the speed of exchange) by changes in the physio- 
logical state, quod erat faciendum. 

Let us now see how justified we are in accepting our results by 
examining the comparison between experiment and independent pre- 
diction. 

1. Since in calculating the k; we have at least on two occasions 
been obliged to use values which were certain only up to an order of 
magnitude, we may reasonably require only that the theoretical and 
empirical decay constants agree up to an order of magnitude. This 
is certainly the case:°® 


Empirical Theoretical 
Ko 96) gL, 
k .058 .03 
Ke ; .02 .003 


We suppose that it would be improbable indeed that ten reasonable 
values assumed for ten independent quantities should yield reason- 
able answers by sheer accident.? 


2. Not only is the theoretical form (exponential sum) the best 
fitting empirical expression, but also the predicted number of expo- 
nentials (n + 1 = 2 + 1= 83) turns out empirically to be the best 
number. 


3. It will appear below that the size, direction, and nature of 
the theoretically deduced partial derivatives is in good accord with 
what one would expect on physiological grounds and/or on the basis 
of experiment. 


The three foregoing observations we have taken to favor the ac- 
ceptance of our equation (1), and hence to warrant the consideration 
of the derivatives in Figure 2 as expressing approximately the effec- 
tiveness of various influences in governing the speed of inert gas ex- 
change. Study of these derivatives permits the following remarks: 


1. The early stages.of absorption are governed chiefly by the 
state of the blood, the intercellular fluid, and the charactertstics of 


6 It has been customary in this field to express decay constants in min-1; we 
have therefore converted from cgs units. a 

7 Where comparison is possible our empirical values are well within the range 
obtained by others. Since the left member of equation (3) is continuous for 
k,; # 7;, we can be sure that normal deviations of physiological values from the 


standards we have assumed will lead to decay constants which themselves remain 
within the normal range. The interesting ordering relationship between the 
k; and 7; has itself important applications (3, 5). 
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aqueous tissue (e.g., muscle), in the sense that these are the influ- 
ences to which ky is most susceptible. 


2. The later stages of absorption are governed chiefly by the 
interaction of the gas with fatty tissue. Thus the smaller the amount 
of fat, the greater the fat solubility, or the greater the fat perme- 
ability, then the more rapidly will the absorption be terminated, for 
these are the influences to which k, is most susceptible. 


3. Both conclusions 1 and 2 can be reached without the untenable 
concept formerly held that each exponential was determined solely 
by the constant of one type of tissue. 


4. Experimental procedures designed to affect one or another 
of the physiological parameters could eventually lead to a rather fine 
control of the absorption curve. 
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CONTRIBUTION TO THE THEORY OF DISCRIMINATION 
LEARNING 


TERRELL ’L. HILL AND LAURA E. HILL 
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In the first two sections, which deal, respectively, with simple and 
double (or successive) discrimination, a comparison is made between the 
theory presented and certain experiments on time discrimination. Sec- 
tion III sets forth a possible theoretical approach to multiple choice dis- 
crimination. 


I 

In a discrimination problem in which there is a choice between 
a “correct” response and an “incorrect” response, we may designate 
the probability of a correct response being made on the n-th trial as 
a(n). In the experiments to be discussed in this section there is only 
one correct response and one incorrect response possible per trial. 
Also, a(0) = 4 because of the symmetry introduced. We are inter- 
ested in the nature of the function a(7). 

We introduce a more or less intuitive postulate whose justifica- 
tion will rest in the experimental verification of deductions made 
from it, rather than in an attempt to deduce it from more funda- 
mental considerations. It may be remarked that the postulate will 
lead to an equation which is a special case of a more general equa- 
tion obtained by H. D. Landahl (1941) from a considerably more de- 
tailed argument. However, the present very simple derivation may 
also be of some interest. 

We postulate that the increase in the probability of a correct 
choice, as a result of additional experience, is proportional to the 
amount of improvement in this probability still possible and to the 
amount of additional experience. That is, 


da=k(1—«)dn, (1) 


where k is a proportionality constant and is a function of the indi- 
vidual’s aptitude and of the intrinsic difficulty of the experiment. In- 
tegrating equation (1) with a(0) = 4, we obtain 


a(n) =1—fe*. (2) 
107 
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Now a is not an experimentally determinable quantity. The nearest 
we can come to it is the average number of correct responses made 
over a certain range in n. However, E(n), the total number of er- 
rors made during n trials, is easily obtainable. We have 


bale, (3) 
dn 
and 
E= f} (i adn (4) 
Hence 
E (n) =~ Bieri wh (5) 


This is a special case of H. D. Landahl’s (1941) equation (34). 


* HALF OF THE RATS WERE NOT USED 
FOR THESE POINTS 


ie 
200 250 300 350 400 


FIGURE 1 


The nature of equation (5) is shown in Figure 1 for different 
values of the parameter k (smooth curves). Obviously, the line 
E = 4 corresponds to k = 0 and no discrimination, while the line 
E = 0 (k= o) holds for perfect discrimination from the outset. 

A. C. Anderson (1932) conducted an extensive experiment on 
time discrimination in white rats in which each rat had to learn 
which of two doors was associated with a short detention period and 
which with a long detention period. Different differences and ratios 
of detention periods were used on different groups of rats in order 
to test Weber’s law. It was found that Weber’s law was indeed satis- 
fied. That is, the average rate of learning was the same, within the 
random error to be expected for relatively small groups of rats, for 
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groups which had the same detention or discrimination ratios but 
different differences. 

Four discrimination ratios were used by A. C. Anderson, r = 4, 
3, 2 and 1.5, with 4, 16, 32 and 8 rats, respectively, as subjects for 
the different ratios. (The discrimination ratio r is defined as 
ry = 1,/t,, where t; and t, are the duration in seconds, say, of the long 
and short detention periods, respectively.) Average values of E were 
calculated from Anderson’s data at intervals of n = 25 and these 
average values are plotted in Figure 1 against n. The smooth curves 
are theoretical curves, from equation (5), for the values of k spe- 
cified in the figure. The points are experimental points. It can be 
seen that the agreement is quite good. 

We now turn to a more detailed examination of the constant k, 
which may be called the discrimination coefficient. It is seen from 
equation (1) that the rate of learning is proportional to k. This 
coefficient may be considered a function of two quantities r and a, 
where r is the discrimination coefficient and a depends on the indi- 
vidual and the difficulty of the experiment (without regard to r). It 
is of considerable interest to discover the form of the function k(a,r), 
since this would make it possible to express Weber’s law quantita- 
tively. 

As we are using average scores in connection with the Ander- 
son experiment, it may be expected that a is nearly the same for 
the different values of r (it would, of course, be advantageous to have 
larger numbers of animals). Hence, plotting k(a,7) from the experi- 
mental data against 7, as in Figure 2, should give approximately, at 


FIGURE 2 ~ 


least, the type of devendence of k on r at constant a. No information 
is available on the variation of k witha. 
Figure 2 suggests* that k(a,r) may have the simple form 


* We are indebted to Professor H. D. Landahl for pointing this out. 
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=a log 7, which agrees with the quantitative formulation of Web- 
er’s law suggested by Fechner. Future and more detailed experiments 
and analysis are necessary to verify this, and to investigate the gen- 
erality of this equation. 


II 

J. T. Cowles and J. L. Finan (1941) performed a time discrimi- 
nation experiment on white rats in which a given rat was first de- 
tained in a chamber for either a short or a long period of time. On 
release, the rat was allowed the choice of two doors. One door was 
associated with the short time interval and one with the long. The 
“correct” response was to attempt to enter the door associated with 
the detention period used on the particular trial. The reward for a 
correct choice was an open door leading to food. Punishment for an 
incorrect choice was a blocked door, which necessitated, in order to 
reach food, a reversal of choice and a passage through the correct 
door. Cowles and Finan point out several faults in Anderson’s experi- 
ment, as an experiment in pure time discrimination, and these faults 
are corrected by the procedure Cowles and Finan adopt. On the other 
hand, although it is true that the discrimination in the Anderson ex- 
periment is a sort of combined space-time discrimination, whatever 
type of discrimination it is, it involves only a single discrimination 
on the part of the animal, while the Cowles and Finan experiment is 
in reality a double discrimination problem. That is, in the Cowles and 
Finan case the rat is first faced with a pure time discrimination, but 
instead of the reward for a successful discrimination being food and 
the punishment an absence of food, the rat receives food in due course 
in either case. Hence, the reward is essentially a short path to food 
while the punishment is essentially a long path to food. Thus, in or- 
der for the rat to appreciate the reward he must be able to discrimi- 
nate between the short and long paths (or time delays). This prob- 
lem involves, then, a double discrimination: a temporal discrimina- 
tion at the outset and then a spatial (or mixed spatial-temporal) dis- 
crimination in connection with the reward-punishment system. The 
object in this section will be to consider such a double discrimination 
problem from a theoretical point of view. 

We shall derive the probability «a that an animal will make the 
correct choice in the following problem: a primary discrimination 
problem is presented and then a choice of two alternatives; the two 
alternatives (reward, punishment) are not of an all-or-none char- 
acter but involve a further discrimination. If the rat performs the 
first discrimination successfully he must also have discriminated pre- 
viously between the reward and punishment in order to choose the 
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correct alternative with surety. As a good example of this sort of 
problem, one might vary with the Cowles and Finan experiment as 
follows. After the initial detention, both doors are left open but on 
passing through either door the rat is detained in a new chamber for 
a certain length of time before being released to food. A long deten- 
tion period is punishment and a short detention period—reward. The 
reward-punishment system here would demand a discrimination of 
the type used by A. C. Anderson and by C. F. Sams and E. C. Tol- 
man (1925). 

Let a; be the probability that a correct choice would be made if 
only the first discrimination were necessary (i.e., with an all-or-none 
reward-punishment system). Let a, be the probability that a correct 
choice would be made if only the second discrimination were neces- 
sary (i.e., if the first “discrimination” were of an all-or-none charac- 
ter). Since a, and a, involve only single discriminations, the equa- 
tions of Section I should hold with 


i 1 = d e-kin é 
ie oar (6) 


In equations (6) 


. k, =k, (a,71), 

and 

Keg = ky (M2 , T2), (7) 
where 7, and 7, are the two discrimination ratios used. In the first 
case mentioned above a, is the probability of the correct choice being 
made but it is not the probability of the first discrimination being 
made, since even when the discrimination is not made the probability - 
of making the correct choice is equal to the probability of making the 
incorcect choice. Therefore, if $e” is the probability of an incorrect 
choice, twice this quantity, or e”, is the probability of the discrimi- 
nation not being made and 1 — e™", the probability of it being made. 
Now, if the first discrimination has been made the probability of 
making the correct choice in the presence of the second discrimina- 
tion is a. , and the probability of not making it is 1 — a,. If the first 
discrimination has not been made, the probability of making the cor- 
rect choice is 4, regardless of a. Therefore, adding the two contri- 
butions toa and1—a, 


a=a,(l1—e™") + be", (8) 


and 
1—o=(1 a) (1-—e*%) +4e™., (9) 
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Hence, 
E=[ (1—a)dn 
wegen) ees (10) 
Dis 2h 
— gS oat a (1 — @- (kat ka) ) | 
2 (key lez) 


We may now return to the experiment of Cowles and Finan. The 
averages given by these authors for six rats at intervals of n = 60 
are used in Figure 3. The dotted curve, which does not seem to fit 
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FIGURE 3 


the experimental results, is based on the assumption of single dis- 
crimination [equation (5)]. This curve was chosen to pass through 
the last experimental point (k = .00086). Other choices would not 
lead to any better general agreement. The solid curve in Figure 8 is 
based on equation (10) and therefore assumes double discrimination. 
Quite arbitrarily, in the absence of other information, it was assumed 
that k, = k,. The curve given is for k, = k, = .00204. Actually, the 
goodness of fit seems to be rather insensitive to the value chosen for 
k,/k2. In an experiment designed intentionally as a double discrimi- 
nation experiment, k, and k, could be obtained independently by let- 
ting 72 > o andr, > o, respectively. Also, if E (cc) is obtained ex- 
perimentally, this provides a relation between k, and k.. 


Ill 
In this section a different type of generalization of equation (2) 
than that considered in Section II will be treated brieflly. As far as 
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We are aware, data are not available at the present time for a com- 
parison of the theory of this section with experiment. 

Suppose that in a single discrimination problem there are N 
choices instead of two. Thus, in an experiment similar to the A. C. 
Anderson experiment, there might be three different detention peri- 
ods, each associated with one of three doors (and detention cham- 
bers) placed in equivalent positions (at the vertices of an equilateral 
triangle). In part of Anderson’s experiment four doors were actual- 
ly used, but they were not placed in equivalent positions. In fact it 
is not possible to place four doors in equivalent positions in two di- 
mensions, since there are three essentially different assignments of 
four different detention periods to the four corners of a square (there 
is only one such assignment for an equilateral triangle and three de- 
tention periods). Hence, for N > 3, experimental results would be 
somewhat questionable. 

Let the detention periods (if the discrimination is time discrimi- 
nation) be t,, &,---, ty in increasing order. Let a; be the probability 


of the door associated with t; being chosen. Clearly, when n = 0, 


(11) 


Oy — Ay — *** = Gy — 
and when n = o, 
o—1, Se = Ag = — Oy —— OF, (12) 


In this general case the correct response is the choice of D, (the 
door associated with t,), and the choice of any other door is incor- 
rect. We postulate here that equation (1) still holds in an appropri- 


ate form: 
do, —k,(1—a,)dn. (13) 


From equations (11) and (13), 


sin) = 1 = ——= Fen (14) 


The discrimination coefficient k, will be a function of t,, t2,-:- , ty. 

Let p; = a, be the probability of D, being chosen, p, be the prob- 
ability of D, being chosen if D, is not chosen, p; be the probability of 
D, being chosen if D, and D, are not chosen, etc. Obviously py = 1 


and ( 
a=—p WT (1—>p,). (15) 


j=l 
We now postulate further that 
dp; =ki(1— pi) dn; LS42N, (16) 
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giving 
Py pomp ecemees La 7 Ne (17) 
N—i+l 

Equation (17) takes the view that if D,, D.,---, and Dj are not 
chosen, the next most “correct” choice is D, and hence pi, by defini- 
tion, should be of the same general form as a, = p.. The coefficient 
k; is the same function of t;, «+: , ty as hk, is of f,,---, tr. 

We define the experimentally measurable quantity S;(m) as the 
total number of choices of D; in n trials. Then 


Si(n) = ff atnyan. (18) 


Substituting equations (15) and (17) into equation (18), one ob- 
tains on integration 


S A Sed (1— e*™) 
— nN — —_— 47m 5 
‘ Nk, 
i Melee oe (1 Shs @- Kat...+ 1-1) 8) 4 > 1 (19) 
N (ky + --- + Kj) 
N—1 


A eee (1 = E- sts) 9) fs 
N(k, +--+» + kj) 

It would be of considerable interest to test these equations ex- 
perimentally, especially in view of the fact that it is by no means cer- 
tain that the set of postulates adopted here are the correct ones. 

The above discussion uses time discrimination as an example but 
is, of course, not restricted to this type of discrimination. 
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In this paper some aspects of the mathematical biophysics of the 
central nervous system, which hitherto have not been treated, are dis- 
cussed. First, a neurobiophysical mechanism for consciousness is sug- 
gested. It provides for the possibility of conscious and unconscious re- 
actions. Next, a mechanism of memory, both on the conscious and sub- 
conscious level, is suggested. The gradual forgetting of remote events 
is ascribed to the inhibition of older memory traces by the more recent 
ones. On the average, an exponential decay of memory with time is 
thus obtained, although memory for unusually strong experiences fol- 
lows a somewhat different law of decay. A homeostatic mechanism is 
then considered which regulates the level of accumulated excitation or 
inhibition. Such a mechanism, under certain disturbing conditions, will 
result in periodical fluctuations of the total cortical excitation with 
periods varying within a very wide range. Finally, a mechanism for 
foresight and desire of future events is suggested. The latter provides 
for the possibility of the formation of subconscious reactions and habits, 
which may be abolished by bringing them into consciousness. 


The recent developments of mathematical biophysics have ac- 
counted for a large number of phenomena in the central nervous sys- 
tem, ranging from a quantitative description of reaction times to 
such complex phenomena as learning (Householder and Landahl, 
1945). An approach to the problem of emotions through the study 
of the mechanism of pleasantness or unpleasantness of visual (Ra- 
shevsky, 1938a, b; 1940; 1942a, b; Rashevsky and Brown, 1944a, b) 
and acoustical (Rashevsky, 1941) stimuli has also led to correct pre- 
diction of new experimental results. The field of the theory of the 
central nervous system is, however, as broad and complicated as the 
manifoldness of the functions of the brain. In spite of having covered 
a very wide range of phenomena in this field, mathematical biophysics 
has still left out much more than it has covered. In the present paper 
we shall outline some approaches to important phenomena in the 
brain, which thus far have not received adequate treatment at the 
hand of the mathematical biophysicists. 

We shall begin by considering neurobiophysical mechanisms 
which may throw some light on the nature of consciousness and mem- 


ory. 
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Let a number of afferent fibers J (Figure 1) lead to the cortex 
and synapse each with a fiber IJ, as well as with the usual cross-in- 
hibitory fibers IJJ. For simplicity, not all of fibers [JI are shown in 
the Figure. The scheme thus far is similar to the one studied before 
(Rashevsky, 1938a, chap. xxii). If the stimulus pattern affecting the 
peripheral ends of fibers J is pleasant, then the same pattern will be 
duplicated in fibers IJ (Rashevsky, 1940, chap. xii). For unpleasant 
stimuli some or all of the fibers IJ will not be excited because an un- 


pleasant stimulus results in a net inhibition of the set of synapses 8, . 
Let, however, each fiber J send an excitatory fiber IV to a different 
cortical level in which it excites a self-circuiting neuron (Household- 
er and Landahl, 1945) Eng, which has a low permanent excitation 
threshold, so that a moderate or even weak intensity of excitation E,y 
of fiber JV throws it into a permanently excited state. Let each Eng 
send neurofibers I’ downwards, so as to produce an excitation of the 
efferent neurons Eff, synapsing with JJ at s,. With this arrange- 
ment even an unpleasant stimulus pattern results in an excitation of 
the efferent and of IJ, as should be the case. Without this or a simi- 
lar arrangement we would be led to the absurd conclusion that an un- 
pleasant stimulus pattern cannot produce any reaction at all. With 
the present mechanism there is, however, an essential difference be- 
tween pleasant and unpleasant stimuli. In the case of the former 
ones, the efferent neurons Eff are activated both through JJ and 
through Eng , while in the latter case the Eff neurons are activated 
- through Eng only. The biophysical correlate of unpleasantness is in- 
hibition of neurons J]. 

Each Eng neuron provides a memory trace or engram of the 
corresponding stimulus, since Eng remains excited after the stimulus 
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subsides. At first glance, one might be inclined to assume that we 
perceive a stimulus pattern directly when the efferents Eff are ex- 
cited through IJ , while we only remember a stimulus pattern if the 
afferent region is excited only through Eng. With this simple pic- 
ture, however, a very unpleasant stimulus pattern which does not 
excite any neurons JJ would not be perceived as real at all, but only 
as a memory. This hardly corresponds to facts in general. However, 
this may have something to do with the feeling of unreality which a 
very painful situation produces in some individuals. 

We now complicate the picture by having excitatory fibers, per- 
haps rather weak ones, connect each neuron J with the correspond- 
ing efferent neuron Eff, by-passing the corresponding neuron JI. 
The total excitation or inhibition of the neurons JI still remain the 
measures of pleasantness or unpleasantness, but even the most pain- 
ful situation will now give a feeling of reality because it will have 
excited the neurons Eff not through Eng only. 

With this mechanism the memory of an unpleasant event is in no 
way different from that of a pleasant event. Neither memory involves 
neurons JJ. This in itself does not necessarily contradict any obser- 
vations. We frequently think without particular pain of past events 
which at the time of their occurrence were extremely unpleasant. 
However, sometimes memories themselves are very painful. To ac- 
count for this we must complicate somewhat our neural mechanism. 


FIGURE 2 


Let each Eng (Figure 2) send “down” a fiber J’, synapsing at 
8’, with a fiber IJ’, which in its turn synapses at s. with the same effer- 
ent neuron as the corresponding neuron JJ. Let each I’ send cross- 
inhibitory fibers III’ to each s’; , so as to form a sort of mirror image 
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of the system of neurons J and IJ. Let, moreover, each neuron J and 
each neuron J’ send excitatory fibers to the corresponding efferents, 
by-passing neurons IJ and IJ’ correspondingly. . 

A directly perceived stimulus pattern is pleasant if it results in 
the excitation of neurons JJ and IJ’. A directly perceived stimulus 
pattern is unpleasant if it results in an inhibition of both neurons 
II and IJ’. A memory of a pleasant stimulus pattern exists when 
neurons JJ’, but not IJ are excited, while a memory of an unpleasant 
pattern exists when neurons II and III are not excited at all, but 
neurons JI’ and III' are excited in such a way that the net result is 
an inhibition at the synapse s’; . 

In general, the memory of a pleasant event will be less pleasant 
than the event itself, since only neurons JJ’, instead of both JI’ and 
II , are excited and excite neurons Eff. The memory of an unpleas- 
sant event is less unpleasant than the event itself. This seems to 
agree with observations. 

Events of which we are conscious we do usually remember to a 
greater or lesser degree. We speak of an unconscious act, if we do 
not ordinarily remember it. This suggests that the excitation of the 
Eng neurons corresponds to consciousness of an event. 

If the synaptic delay along fibers IV , each of which may actually 
be a chain of fibers, is long, then some stimulus patterns may pro- 
duce a reaction through the neurons Eff, before the neurons Eng 
are stimulated. Such a reaction is subconscious and is not consciously 
remembered. If each of the efferent neurons Eff sends an inhibitory fi- 
ber V to one of the synapses of the corresponding chain of fibers IV, 
then once such a subconscious reaction is produced, before the corre- 
sponding neuron E'ng is excited, a subsequent excitement of that neu- 
ron Eng will be prevented. The act remains subconscious. In the 
absence of such an additional mechanism, or when the reaction is 
short, and therefore the inhibition of the chains JV is short, the neu- 
‘rons E'ng may be stimulated after the corresponding efferent neurons. 
Thus we may or may not remember afterwards an act of which we 
were not conscious at the time of its performance. 

It is, however, known that subconscious acts may in a way be 
remembered also, though the memory itself is on the subconscious 
level. Under the action of some drugs, such as sodium amital, such 
“subconscious” memories may be evoked. The above picture does not 
account for that, and we must complicate it further. 

Let each of the fibers V be a chain of fibers, and let that chain at 
some place s;y synapse with a self-circulated neuron Eng, . Let each of 
the chains of fibers V, mentioned above, inhibit through a last inhi- 
bitory link a synapse of the chain IV somewhere between Eng, and 
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Eng. Let, moreover, the chain of fibers V have at one of its synapses 
a self-circuited circuit C, whose permanent excitation threshold is 
such that when V is excited weakly the excitation is reversible, but 
when V is excited strongly, the circuit is excited permanently and 
the part of the chain V between the circuit and sy (Figure 2) remains 
permanently excited. 

With this mechanism, once a sufficiently strong subconscious re- 
action is produced, the corresponding Eng, are excited but not the 
corresponding neurons Eng. This may be interpreted as meaning 
that the individual will not consciously remember the event but it 
will be stored in the subconscious memory. If through the action of 
some factor, for example sodium amital, the chain V between C and 
Sy is inhibited, the excitation of Eng, causes an excitation of Eng, 
and the event is remembered. 

The individual’s reaction to such a “disinhibited” previously sub- 
conscious memory will be through Eng, I’, s’,, and II’. In this case, 
once a subsconscious memory is evoked into consciousness, it remains 
conscious. The chain JV may, however, branch off at sy; an excitation 
chain VJ which connects directly to Eff. In that case, disinhibition 
of sy may produce a reaction through VI and Hff, without neces- 
sarily exciting the neurons Eng. The reaction may be interpreted as 
the individual’s statement about a past event remembered only sub- 
consciously. But when the disinhibition is over and the neurons Eng 
did not become excited, the individual will again not remember the 
event. All will depend on the intensity of excitation of Hng, or of 
the intensity of excitation H,, of the neuron chain between Eng, 
and Eng. If Ep is greater than the threshold h* of the permanent 
excitation of Eng , the inhibition of V will result in a conscious mem- 
ory of the event, even after the disinhibition is over. Otherwise, the 
individual will again forget consciously the event. 

It is empirically known that memory of an event becomes weak- 
er as time goes on. One might attempt to interpret that in terms of 
a gradual extinction of the excitation of the engramic neuron. Such 
an interpretation meets, however, with a fundamental biophysical dif- 
ficulty. The time factors for all kinds of processes in single neurons 
are of the order of small fractions of a second. Memories, on the 
other hand, may last for years. Hence we would have to introduce 
ad hoc some very slow process as an additional hypothesis. Invok- 
ing in this case the possibility of very slow somatic changes does not 
seem plausible. The following scheme takes care of the situation in 
a much more natural way, without any additional assumptions be- 
yond those usually made in the mathematical biophysics of the cen- 


tral nervous system. ‘ 
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Let each engramic neuron Eng send to every other one inhibi- 
tory fibers. In the stream of experience more and more engrams be- 
come excited. Hence each of the already excited engrams becomes 
more and more inhibited by the others. The memory for any event 
will gradually weaken because of inhibition by memories of new 
events. Thus a general weakening of memory will occur, as time goes 
on. But this weakening, contrary to facts, would affect equally the 
recent and the long past memories. Only the strongest memories will 
eventually remain, the weaker ones being completely inhibited. 

The network shown in Figure 2 provides for such a mechanism, 
though in a too restricted sense. According to Figure 2, we should 
not remember too complex situations, or remember several things at 
one time, for all involved engrams will become inhibited when their 
number exceeds a critical value n* (Rashevsky, 1938a, chap. xxii). 
With such a scheme we would find that if on the average a constant 
number of engrams are formed per unit time in the stream of experi- 
ence, then after a finite time all engrams will become inhibited and 
no memory will be possible. 

A different mechanism is provided by cross-inhibitory fibers go- 
ing from synapse to synapse (Rashevsky, 1938a, chap. xxii, Figure 
56), as shown on Figure 3. 


Let E' be the excitation of each engram. For simplicity we now 
assume E' to be the same for all of them. In the stationary state the 
value of ¢ produced at the synapse is, using previous notations (Ra- 
shevsky, 1938a), 


A 
Pere nt rt os (1) 


Denoting by h the threshold of the cross-inhibitory fibers, we find 
that the excitation H; of each of the inhibitory fibers is 
E;=B(e—j—h), (2) 


where j is the total amount of inhibitory factor at the synapse, while 
the amount of inhibitory factor produced by each inhibitory fiber is 


ied 
v] SS i (3) 
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Equations (2) and (3) give 


7 =B yi(e—j—h). (4) 
If the total number of excited engrams is », then at each synapse 
j= (n—1) 7 =B8 yr (n—1) (e—G—h). (5) 
Putting 
a j — ri) , (6) 
we have from equations (1) and (5) 
$=7E—By(n—1)(¢—h), (7) 
which gives 
vE By (n—1) 
o= it ky ere (3) 


—— 14+By(m—1)) 14+8y,(n—1) 


For any finite the excitation ¢ is positive, and tends to h as 1 be- 
comes infinite. If the threshold h is very small, we may neglect the 
second term of the right side of equation (8) and obtain 


(pened 
1+Byi(n—1)° 
Since 7 varies linearly with time, therefore ¢ tends to zero as 

Vt. 

If one engram is excited with the intensity E, and the other 

(x — 1) engrams with intensities ZH, < EH,, then we have, neglecting 


(9) 


, 


&—vyE,, &g= y Ee; 


; : (10) | 
ji=Byi(n—1)d2; fe=Pyi(t—2)bd2ot+ Burgi. 
Hence with 
Byu=n, (11) 
we have 
od. —y Hi — n(n —1) d23 
(12) 
oo = y 2 — 9 (N—2) $2 —7 fr - 
This gives 
peas Heactira St Ba) Ei say 1 Ls ; (13) 
7 71 a) 
E.— E, 
bo= i vn (14) 


(l= 9)? +L —9) 0 
For FE, = E, those equations reduce to expression (9). 
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If E, and E. are of the same order of magnitude, then for ¢' 
and ¢2 to be positive we must have 7 <<1. If 7 << 1 and E,/EHs > 
1/y, then ¢. < 0. Physically this means that none of the inhibitory 
fibers converging on the engram with excitation EH, are excited. Hence 
in that case ¢, = a H,. Again only the strongest engram is excited 
at all, the others being inhibited. 

The above mechanism results in a weakening of memories with 
time, the weakening being the same for recent and long past events. 
Actually, recent events are remembered better than old ones, other 
conditions being equal. This suggests that while newly formed en- 
grams inhibit those formed previously, they themselves are not in- 
hibited, at least appreciably, by those previous ones. 

The following mechanism provides for such a situation. Let any 
two engrams be connected as represented on Figure 4. Let the in- 
hibitory 12’ and 21’ be of the same strength as the excitatcry fibers 
12, and 21. If J is excited first, then 12 and 12’ will cancel each other 


FIGURE 4 


out and JT will not be inhibited. But if after that JJ is excited, then 
I will be inhibited via 21’ because 21 is inhibited by the previous ex- 
citation of 121’. The constants of 121 must be such that 121 is ex- 
cited strongly even by weak stimuli, so that when the excitation of I 
is reduced through inhibition by IJ, still 121 would inhibit 21. Due 
to the symmetry of the arrangement, if IJ is excited before I, IJ will 
be affected by J, but not I by II. 

Consider now n engrams all excited at different times. Denote 
by 1 the engram which is the last to be excited; by 2 the one before 
last, and so on. Let their intrinsic intensities of excitation be ail the 
same and equal to EZ’. Then for the first engram, we have 


di—ak, (15) 


since it is not inhibited by any others. For the second we have, neg- 
lecting the threshold of the inhibitory fibers, 


do =aHk—nd=—aEk(1—7). (16) 
For the third we have 
$3 =a BH —74(d, + $2) =a (1—2y + y*) =aE(1—7n)?. (17) 
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If 
on — (1 —7)*™, (18) 
then 
dns — a BE — (di + de --+ dn) 

=a #{i—nf{l + (L—y) + (l—n)? +--+ (L—n)™]} 
(19) 
i (l—7n)*—1 

| 


Since equation (18) holds for n = 1 and 2, it holds for any n. 
If x is proportional to the time ¢ , then we have 


o:=A(1—n)*, (20) 


which, since 0 < 7» < 0, is a decreasing function of t. Thus the long- 
er ago the engram was formed, the weaker it is. If at some time an 
engram is formed with an intensity E* > E’, and followed by n — 1 
engrams of intensity EH, then 


gn — a E* — ($i + 2 +--+ + haa) = 
af°—akFynfit (l1—7n) +---+ (1—y)™*] = (21) 
ie Oder Dg) OE] 


=akF;1+ =aH(1—y)”. 


As n increases and as time goes on, ¢, decreases but tends to a con- 
stant value a(#* — EF). Thus a few strong memories outstanding in 
the stream of experience will slightly decrease, but then remain con- 
stant. Of course all engrams are of different intensities, but we may 
consider them as of an average same intensity EL, except for a few 
outstanding memories. If on the average outstanding memories will 
occur ever so often, then their intensity will also decrease with time 
to zero, since as time goes on their number mounts to infinity. 

From the discrete case we may pass to the continuous and write, 
counting t backwards, 


el 
o(t) =a H(t) —n | g(é)dé. (22) 
6 
Differentiating, and denoting 
dE 
—=E', 23 
dt (23) 
we find 
do 
——akE'(t) — 5 24) 
wus (t)—n¢ ( 


The general solution of equation (24) is 
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g=et (boty ‘i E' (t) ev dt}. (25) 


If E(t) = const., then E’(t) = 0 and equation (25) gives 
p= poe", (26) 


where ¢» is the intensity of memory just formed. 

We shall assume here that the above discussed mechanism of 
memory extinction holds only for the “conscious” memory. We shall 
consider that the neurons Eng, do not cross-inhibit each other, and 
that therefore subconsciously even the remotest events may be re- 
membered. The number of events stored in the subconscious is lim- 
ited only by the total number of neurons Eng, available. 

In the stream of experience there is a constant increase of exci- 
tation on the level of neurons Eng,. There is thus an accumulation 
of the level of the total net excitation ¢ = « — j in the cortex. It is 
natural to assume some homeostatic mechanism which will tend to 
maintain a constant level ¢,. One of the simplest mechanisms of that 
kind would be the following (Figure 5). The impulses coming from 


FIGURE 5 


the external world are relayed to the cortex C through the thalamus 
T. While we may consider that in an ordinary stream of experience 
not upset by any too unusual events, the average flow of impulses to 
the thalamus will be constant, yet even then the intensity S of the 
individual impulses will be scattered around the average according 
to some distribution function E(S). Let all the Eng, neurons of C 
send fibers to the fibers CT (Figure 5), which are a pair of fibers, 
one excitatory, the other inhibitory. Since ordinarily the changes in 
the flow of outside stimuli are slow compared with the time constants 
of the neurons, we may consider the fibers CT as always in a quasi- 
stationary state of excitation. 

The constants of the two fibers may be chosen so that for weak 
stimuli at C the excitatory fiber has almost reached its saturation 
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value, while the inhibitory one is only weakly excited. For strong 
stimuli the excitation of the inhibitory fiber catches up with that of 
the excitatory, and eventually exceeds it. 

But then for high values of the total amount of 6, = 4 — do ac- 
cumulated over a long period of time, the pair CT will act as an in- 
hibitory fiber, raising the thresholds of the fibers TC , running from 
the thalamus to the cortex, and thus cutting off from the cortex some 
of the weaker stimuli of the stream of experience. This will result 
in a decrease of the rate of accumulation d¢,/dt. If the average rate 
of accumulation of ¢,, that is, the flow of stimuli is constant and 
equal to S,, then we have, denoting by a and {4 two constants of pro- 
portionality, 


eee EP 
= a ed 1 . 
ta fe (27) 
By differentiation we obtain from equation (27) 
aby 
de - A $1, (28) 


the solution of which represents undamped harmonic oscillations with- 
in a period 
meee eae 
Va 

In this simple case we would thus have periodical fluctuations of 
cortical excitation and inhibition. Unless at some time a particularly 
high value of S, has brought ¢, to a very high value, the amplitude 
of the oscillations of ¢, will be usually small. But if an external very 
strong experience has made ¢, very large, the amplitude may become 
very large, and abnormally strong fluctuations of the state of excita- 
tion of the cortex will result, such as are observed in manic-depres- 
sive psychoses. 

Let us consider the somewhat more general case when S, is not 
constant but is itself a function S,(t) of time. Then instead of equa- 
tion (28) we shali have, putting dS,/dt = S,'(t), 


dh, 
dt? 

Let in general S, be constant, except for a short time when it rapidly 
raises to some value S,, , then after an interval of time T, again drops 


to the original constant value. The quantity S,’ then will raise from 
zero to a maximum value, drop again to zero, become negative, and 


ik (29) 


+2 ¢,=aS,'(t). (30) 
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again return to zero. It will be represented by a curve as that shown 
in Figure 6. But then for that interval of time, equation (30) repre- 
sents a forced harmonic oscillation. If 7, is very close to 2x/VA, 
then even one “wave” of the external force, S,’(t), may increase the 


FIGURE 6 


amplitude of ¢, very appreciably. An equal temporary increase of 
S, to the value S,, but for a period 7; which is much different from 
2n/\/4 will not result in such abnormally high amplitudes. It must 
be noted that bringing by some external means ¢, and even S, to zero 
at some moment does not necessarily result in the cessation of the 
fluctuations because d¢,/dt may remain at that moment either strong- 


t 
ly positive or negative, depending on the value of the integral if ¢idt . 


Equation (27) should of course be considered only as a first ap- 
proximation. The relation between d¢,/dt and the cumulative value 
of ¢, will in general be non-linear. The additional of non-linear terms 
will introduce harmonics into the oscillations and the possibility of 
periodical very sudden changes in ¢, will be present. 

If we consider another more general case, namely that the Eng, 
neurons inhibit each other, or that the Hng neurons also connect with 
CT , then in the cumulative effect the amount of ¢, which corresponds 
to a more remote time, contributes less than a recent ¢,. If we as- 
sume an exponential decay, as in equation (26), then instead of equa- 
tion (27) we obtain, for a constant S, , 


t 
= aS,—2 { di (r) e-7'* dr, (31) 
Bringing e-”' outside of the sign of the integral, multiplying by e7*, 
differentiating with respect to t, and then cancelling e-7' , we obtain 


Ve dd, 


| de" dt 
which represents damped harmonic oscillations, the amplitude de- 
creasing as e-?’, The greater 7, in other words, the weaker the mem- 
ory, the more rapidly the oscillations will die down. The value of ¢, 
oscillates now not around ¢, = 0, (¢ = ¢o), but around ¢,/=a7 S/d. 


+idi=anSo, (32) 
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The neurobiophysical mechanisms discussed above provide for 
aftereffects of past experiences through memory. We now shall brief- 
ly discuss mechanisms which will have to be studied for the interpre- 
tation of the effects of anticipation of future events upon the indi- 
vidual. Fundamentally, any idea of a future event presupposes that 
the different elements of that event have been experienced in the 
past. The combinations of those elements may be entirely new, and 
never experienced before. But the elements themselves must have 
been given by previous experience. 

It is natural to consider a generalization of the scheme of Figure 
2, in which case each neuron Eng excites a fiber (Figure 7, VII) lead- 


ENGs 


FIGURE 7 


ing to a neuron F in some different region of the brain. While the 
excitation of an Hng neuron corresponds to a memory of a past event, 
a similar excitation of a neuron F' corresponds to the idea of the event 
happening in the future. Since according to that scheme, a present 
event also excites the corresponding F' neuron, we may use the fol- 
lowing interpretation: 

An event is perceived directly when the corresponding peripheral 
afferents, the Eng and F, are excited. If only Eng and F are excited, 
we have the memory of the event. If only F is excited, we think of 
the event in the future. - 

If in the past an experience A has been followed several times by 
the experience B , then the temporal pattern A — B will be “remem- 
bered” through the corresponding neurons Eng. During the experi- 
encing of the temporal pattern A — B, the corresponding neurons 
(or perhaps better, groups of neurons) F', and F's will also be excited. 
If the different neurons F are connected through “conditioning” path- 
ways, discussed elsewhere, (Rashevsky, 1938a, chap. xxv ; Householder 
and Landahl, 1945, chap. xi), then a simultaneous or closely succes- 
sive excitation of F, and F; will have as a result that excitation of 
F, alone will, through the conditioning pathway, elicit the excitation 


of Fz. 
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Let now experience A happen alone. Through the corresponding 
neuron Eng it will produce the excitation of F,. Through previous 
conditioning F’; will become excited. But there will be no excitation 
of either the peripheral afferents nor will there necessarily be a very 
strong excitation of the Eng corresponding to the event B. Hence, 
the occurrence of A will make us anticipate B in the future. We have 
here an elementary mechanism of foresight. 

Such phenomena as rational problem solving or planning are 
probably closely connected with the functioning of neurons F’, involv- 
ing mechanisms outlined elsewhere (Rashevsky, 1938a, chap. xxix). 

When an external event evokes our concern or thought of a fu- 
ture event, it is frequently in the biological interest of the individual 
to concentrate on that future event for a sufficient time. A mechanism 
for such a concentration may be provided by the scheme shown in 
Figure 7. A neuron F excites some other center C’, which in its turn 
excites F’. This enhances the excitation of F. With properly chosen 
constants of the neurons involved, the circuit FC’ will not remain ex- 
cited if the original excitation of F is removed. But if in a pathologi- 
cal individual the constants of the circuit are different, then even the 
removal of the original stimulus A , which resulted in the excitation 
of F, will leave the circuit FsC’, and hence F;, excited. The indi- 
vidual will not be able to bring to a stop his concern about the future 
event F';. If, as is natural to assume in accordance with the general 
mathematical biophysics of the central nervous system, each F' sends 
inhibitory fibers to other parts of the brain, such a strong excitation 
of one group of F' neurons will result in an inhibition of other func- 
tions of the central nervous system. Such a general inhibition may be 
removed by cutting the circuit F';C’. 

It is plausible to assume that the F' neurons are located in the 
frontal lobes. It may be suggested that what is achieved in prefrontal 
lobotomy is just the severing of the circuits F;C’. Once those circuits 
are destroyed, not only does the permanent preoccupation with one 
idea cease, but it becomes impossible for an individual to properly 
concentrate on future events. Complaints of lobotomised patients 
about a certain “haziness” of the future have been occasionally re- 
corded in the literature. 

Let a pattern A + B be several times followed by a stimulus C, 
which produces the reaction R,. Let, however, A alone not be fol- 
lowed by C,. As a part of the pattern A'+ B, A will continually pro- 
duce reaction R,. Let A, presented alone, be followed by a stimulus 
C, which unconditionally produces the reaction R;. While A is ex- 
perienced alone several times, it not only becomes conditioned to R, , 
but (Rashevsky, 1938a, chap. xxv) through it inhibits its own con- 
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nection with R,. With properly chosen constants of the fibers, the 
process of extinction of the conditioned reflex A > R, is slow. More- 
over, after a period of rest when A was not experienced at all, an 
experiencing of A will again result in the reaction R, until A has been 
repeated alone a sufficient number of times. Hence in general, an ex- 
perience of A will produce both R, and R.. If the efferent paths of 
f, and R, cross-inhibit each other, only one of the reactions will be 
produced, namely the stronger one. 

If B has been always followed by C, while A alone has been al- 
ways followed by C,, then we may say that reaction R, is the “cor- 
rect” reaction to B , while the reaction R, is the “correct” reaction to 
A alone. The experience A may have occurred several times in con- 
junction with B accidentally, and thus become conditioned to R,. If 
A — R, has been very strong, A > R, will not overcome it, and A 
will always result in the wrong reaction. If, in the pattern A +B, 
the central neurons, which form the beginnings of the conditioned 
circuits, cross-inhibit each other, then the conditioning of A to R, 
when in conjunction with B is weaker than if A were alone; for that 
reason the reaction A — R, will be weaker than A — R,. Hence when 
A is presented alone, it will elicit the correct reaction R.. 

But if the cross-inhibition between A and B is weak, the reflex 
A — R, may become so strong as to always prevail over A > R,. 

Let A stand for the experience by a person of seeing another in- 
dividual, and let B stand for this individual’s performance of some 
hostile act. R, may then be a fright reaction, or the feeling that the 
other individual hates the person. On the contrary, R, may be a re- 
action of friendliness. If we always have A — R,, then the person 
has the feeling that every other individual is hostile to him or her. 
As follows from the above, this type of paranoia may be caused by too 
weak cross-inhibitory fibers between the afferent neurons. 

The same situation may, however, be caused by an entirely dif- 
erent mechanism. Let R, be a very unpleasant reaction which results 
in a general cortical inhibition. Then if we first have A + B > R, 
with a resulting A — R,, the process of conditioning of A alone to 
R, is inhibited because before the conditioned reaction A — R, is es- 
tablished, the already existing reaction A — R, produces inhibitions 
of all centers. (ee 

If, in particular, we assume that any unpleasant situation results 
in an excitation of a definite “pain” center which sends off inhibitory 
fibers to other parts of the brain, then the above discussed situation 
may be remedied by cutting the fibers leading from the rest of the 
cortex to the pain center. A surgical treatment of some paranoid con- 


ditions thus appears as a possibility. | 
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As we have seen, qualitatively the same type of paranoid symp- 
toms may be due to entirely different mechanisms. The next problem 
is to develop the different quantitative aspects of the two mechanisms 
by means already used in mathematical biophysics. In this way 
through quantitative tests it may be possible to establish which of 
the possible mechanisms is operating in a given case. This would 
have both diagnostic value as well as give clues to possible therapeutic 
procedure. 

Let it be possible to condition Eng, to any reaction, and let the 
Eng’s send inhibitory fibers to all conditioned pathways going from 
the corresponding Eng,’s. Then the mere bringing into consciousness 
of any subconsciously remembered event may abolish some subcon- 
sciously formed undesirable reaction. 

Let now F; be excited by a center F’, which also is connected 
with both the corresponding Eng and Eng, by tracts VIII, and VIII, . 
The excitation of F';’ may be interpreted as corresponding to the de- 
sire of B. Due to the type of connections shown in Figure 7, when 
we desire B , we also zpso facto think of B in the future. But we may 
think of B in the future (excitation of F;) without desiring it. 

Suppose at some time the desire for B was followed by an event 
A , which produces reaction R. Since F;’ through VJI/, and VIII, ex- 
cites both Eng and Eng, , either of which can be conditioned to pro- 
duce R, R will be produced later on whenever the desire for B is 
experienced. If the excitation Evm, of VIII, is always stronger than 


Evin, of VIII, ; then the “subconscious” part of the conditioning, that 


is, the conditioning of Eng, to B, will be inhibited by the excitation 
of the corresponding Eng. Thus the reaction R will be produced 
quite consciously, and if it is an undesirable reaction, it may be con- 
sciously suppressed. But if in an individual the connection VIII, is 
always much weaker than VIII, , then the reaction R will be produced 
subconsciously. Reaction R may be undesirable, even antisocial, or 
even painful to the individual. Thus B may be attention by others, 
while A may be some physiological disturbance or mild illness pro- 
ducing the reaction R, say a headache or other symptom. If desire 
for attention (B) has been accidently followed immediately by A , and 
thus produced R , R may be produced again every time B is desired, © 
although A (the actual physiological disturbance) may be absent. We 
have here a mechanism that produces some simple forms of hysteria. 
A hysterical individual, according to that picture, is characterized 
by Evin, <E vin, + 


From this point of view, severing of VIII, may abolish some 
types of psychoneurotic behavior. In this connection it may be inter- 
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esting to speculate concerning the mechanism of the beneficial effects 
of prefrontal lobotomy. If all connections to and from the frontal 
lobes are severed, this would mean severing VIII , as well as the cir- 
cuit F; — C’ (Figure 7). The severance of the latter would, as men- 
tioned above, abolish some types of depressions connected with con- 
stant concern over the future. Severance of VIII, would abolish cer- 
tain types of hysterical behavior and subconscious habits. It may be 
worth noticing that improvements in psychoneurotics have sometimes 
been reported after lobotomy. Perhaps with further experimental 
study and improvement of surgical technique, it may be possible to 
sever at will either V/II, or circuit Fs — C’. Different kinds of “par- 
tial” lobotomies may perhaps be used in the future for different types 
of mental disorders. That a complete severance of all connections may 
not be necessary, seems to be indicated by some recent work (Hof- 
statter et al., 1945). 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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Developing further a recent suggestion of H. D. Landahl, equations 
are derived which give the probability for an individual to judge two 
stimulus patterns as similar or dissimilar. The possibility of experi- 
mental verification of those equations is discussed. Next, a mechanism 
is described which provided for abstraction by responding only to “es- 
sential’ components of a stimulus pattern. Equations, verifiable in prin- 
ciple experimentally, are derived. Finally, a mechanism is suggested 
for logical inferences, and equations are derived which give the prob- 
ability of making an error in a reasoning consisting of a chain of syl- 
logisms, as well as the probability of being unable to complete the chain 
of reasoning at all. 


Recently H. D. Landah] (1945) has suggested a neurobiophysical 
mechanism which provides for the reaction of judgment of similar- 
ity and dissimilarity of two multimodal stimulus patterns, or in com- 
mon parlance, of two objects. He suggested some applications to 
pathological cases of erroneous identification and confusion. The pur- 
pose of the present paper is to develop H. D. Landahl’s suggestion 
further and to discuss a similar mechanism which provides for ab- 
straction and for some forms of logical thinking. We presuppose the 
reader’s familiarity with H. D. Landahl’s paper. 
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We first shall somewhat simplify H. D. Landahl’s mechanism for 
perception of similarity. Let a mechanism A (Figure 1) for percep- 
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tion of differences in one modality, such as suggested by H. D Lan- 
dahl, inhibit through the fiber VIII’ a synapse of a chain of fibers 
VIII, which produce a reaction R,, the similarity or identity reac- 
tion. The chain VIII may either be excited by a permanent self-ex- 
citing center or, what is more natural, may be excited by either of the 
two stimuli acting on the “difference mechanism”. In this case, if the 
difference A; = 0, then no “difference reaction” Ra is produced, but 
the similarity reaction R, is produced. If A; > 0, then Ra is produced, 
but R, is absent. 

According to this mechanism, similarity is a purely negative con- 
cept meaning absence of difference. A situation, mentioned by H. D. 
Landahl, that two things may be neither similar nor different, may 
occur with the present mechanism also, if the chain of fibers VIJIJ is 
broken or inhibited somewhere between the synapse s; and the final 
path. Such a situation would, however, be pathological, whereas with 
Landahl’s mechanism it may occur normally. 

We now shall compute the probability p judgment of similarity 
in terms of 4;. The probability q of a judgment of difference is sim- 
ply 

C= Lame (1) 


It follows from H. D. Landahl’s discussion (1945) that the prob- 
ability » is the same as the probability of a “doubtful judgment” in 
the mechanism for psychophysical discrimination studied by H. D. 
Landahl previously (1938). We may take over directly some of the 
equations developed by Landahl (1938). Since our A; is essentially 
the same as ¢,; — e, of Landahl (1938), we have, assuming a distribu- 
tion function p(&) for the fluctuation of excitation of the form 


p(é) Fen, (2) 
the following expressions (Landahl, 1938, p. 116): 
for A= A: 

p=—1—e"coshkA, (3) 
q=e*"coshk A. (4) 

For d= h: 
p=e(lL— ey. ~o = 21 + om), (5) 

For A> h: 
p = e*4 sinh kh, (6) 


q=1—e*sinh kh. (7) 
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Equations (3) and (4) as well as equations (6) and (7) reduce to 
equations (5) forA=h. 

The use of expression (2), while having the advantage of giving 
simple closed final expressions, has the disadvantage of requiring the 
discrimination between the three cases 4 << h, A=h,andA>h. A 
convenient substitute for expression (2), which approximates a nor- 
mal distribution rather well, especially for small values of &, is* 


2 


p(é) =A(#—ae+—), (8) 


“= 


used for values of é in the interval — V/a/2, + \/a/2, A and a being 
constants. It is readily seen that within the above mentioned interval 


FIGURE 2 


the curve p(é) has the shape as shown on Figure 2. The use of ex- 
pression (8) amounts to neglecting the very small probabilities, given 
by a normal distribution function, for large values of €. Instead of 
the requirement 


[ p@de=1 , (9) 
we now shall require 
+ a/2 at 
af sot ees 1 (10) 
-ya/2 
Requirement (10) leads to 
A=1 8407", (11) 
so that finally we have 
p(é) = 1.84 a9? (£4 —a & Dale (12) 


The probability p of similarity reaction is given by (Landahl, 
1938; Rashevsky, 1940) 


* Suggested by L. L. Thurstone. Personal communication to the author. 
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p= [2 dé. (13) 


(A+h) 
If expression (12) is used, we must always have 
AGE WN aie, ASh or A—if., (14) 


With this restriction, we find, by introducing expression (12) into 
equation (13) and integrating, 


2 2a 
p= 1.544°° E At — 2(a — 2h?) A’? ja PSE eae’ a1 a) ee NTE) 
'-) |\22 Se 


G, =A By 
FIGURE 3 


Since, because of inequality (14), a — 2h? > 0, therefore p varies as 
a function of A as shown on Figure 3. By solving the equation 
0p 
—=0, 16 
ae (16) 


we find that the curve has a horizontal tangent for 


4 aes ie - 
2 2 


Equation (15) has a meaning only in the interval (0, 4,) since for 
A > A,, », as given by equation (15), increases while it actually 
should drop monotonically to zero. For A = A,, p is positive, but 
very small. 

For later purposes we shall now derive an expression for h in 
terms of the thresholds of the different fibers involved. 

If, with H. D. Landahl (1938), we assume that the excitatory 
and inhibitory fibers have the same parameters, that is, A—=B,a=b, 
then the intensities E,;, and E,,’ of the fibers I/I on one hand and IJI' 
and IV’ (Figure 1) on the other are given by 


En = Be: — hur) 5 Ey’ = B(e2— hr’). (18) 
Hence, denoting by «, and j,, the values of « and j at the synapse 8, , 
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AB erp 
SS Se (e, — hur) ; Ja = —_ (& — hy’). (19) 
a a 
Therefore 
| ak ABs 
&3 — Ja — — (8; — &) + — (hy — Ain) , (20) 
a a 
and similarly, at the synapse s, 
| Ap 
&, — Js = — — (8,—&) + — (Ayy’ — Rin) « (21) 
a a 
We do not have in general 
&3 = ja —— (8,— 3), 


as in Landahl’s case, unless hyy’ = hy, , that is, unless the thresholds 
of the excitatory and inhibitory fibers are the same. In the general 
case it may happen that both e; — 7, and e, — 7; are positive and great- 
er than hy. This would complicate matters considerably. However, 
if 


AB 
0< Sat (Ary’ al hin) = hy, (22) 


then, with «, — e. > 0, we always will have 
fee Us &—js<h,. (23) 


Assuming, therefore, inequality (22) to hold, we now have for the 
condition of a response Ra, 


é,—ji>h, 


which, together with expression (20), gives 


A : 
i oo 3B hyt+thu—lhy. | (24) 
Hence 
a ; 
TiGian Keukh os : | (25) 


The greater the threshold of the cross-inhibitory fibers, the smaller 
h. The quantity h also increases with the threshold hy . 

Now we proceed with computing the probability of a similarity 
or dissimilarity judgment for two multimodal ppaee patterns. We 
may consider two different cases. 
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Case A. For each modality there is a separate neuron chain for 
a reaction R,. The whole mechanism consists of a number of such 
mechanisms, as represented on Figure 1. Let the subscript 7 refer to 
the t-th modality. Then if all Ra,’s are zero, all R,; are produced. If 
only some Ra, ’s are zero, then only some of the R,; are produced. With 
this mechanism an R, reaction is produced whenever at least one Ra, 
is zero. In other words, an individual will always notice a similarity 
whenever there is one in any respect. 

The probability of a given A; to give no difference response is 
9(A;). The probability of a given A; to give no difference response, 
but for all other 4;’s to give one, is 


p(A;) q(Ar); (k2%). (26) 
k 


Hence the probability that any one A; will give no difference response 
while all others will is 


Sp(As) Wq(4); (kt). (27) 
$=1 k 
The probability that a given pair A; and A, will simultaneously 


give no probability response, regardless of what the remaining A,’s 
do, is 


p(Ai)p(A:). (28) 


The probability that the above mentioned pair gives a no difference 
response while all other A;’s do give a difference response is 


p(Ai)p (Ar) IT q(Ax) ; (kK2i;k 21). (29) 
k 


Hence the probability that any one pair gives a no difference 
response, while all other A;’s give a difference response, is 


2X DP(A:)p(Ar) Wal); (kK2t;k21;121). (30) 
| ae k 


Similar expressions are obtained for the probability of any three 
A’s giving a no difference response, all others giving a difference re- 
sponse, etc. Summing all those probabilities, we finally obtain for 
the total probability of a similarity response 


P= (Ai) atghcia) +> 2 P(Ai)p (Ai) I (Ax) 
a fist S k 


+o BoE (Ay) p(4,) Hl g(Ay). er 
i r k 
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Case B. There is one common response path for R, for all mo- 
dalities (Figure 4). Now R, is produced only when all Ra ’s are 
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FIGURE 4 


zero. Hence the total probability P of R, is equal to the probability 
of all A;’s giving simultaneously no reaction Ra. This is equal to 


P=II p(A;). (32) 
t=1 

Using expressions (3), (5), and (6), or expression (15), for 
p(A;), we can compute P as a function of all A,’s. 

In order to compare equations (31) and (32) with experimental 
data, it is necessary to have a series of stimulus patterns with known 
values of A;’s for each pair. For some modalities the differences 4,’s 
may be measured directly by the difference of some physical quan- 
tities. Thus we may have a set of objects of different size, weight, 
temperature, and perhaps saturation value of a given color, or still 
better, intensity of illumination. In some cases, however, the modal- 
ities may not be of a simple physical nature. Thus we may compare 
two objects in regard to the degree of their difference in being used 
to perform a certain function. For instance, a bell and a horn would be 
generally considered as more similar in regard to their use as sound- 
ing devices than a bell and a dishpan, although the latter may be used 
for sound-producing purposes. In such cases a scale of A; values 
may be established by standard psychophysical methods, for instance, 
the rank order method or the method of paired comparisons. We may, 
for example, ask a large number of subjects to arrange a number of 
objects first in order of their use for a certain function, then in order 
of values of some other modality, etc. From such rankings scale val- 
ues for different modalities may be computed. Then pairs of those 
objects may be presented to one test subject many times, with the 
request to pass a judgment of similarity or dissimilarity in general, 
and from the observed probabilities of different responses equations 
(31) or (32) may be compared with experiment. 

Consider a stimulus pattern which becomes conditioned to a def- 
inite motor reaction, in particular to a verbal reaction. The verbal 
reaction is then the “name” of the stimulus pattern, and stands as a 
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symbol for the latter. 

In general, stimulus patterns, denoted by the same word or hav- 
ing the same name, are somewhat different with respect to some of 
the component single stimuli, but they all have some component stim- 
uli in common. Thus a “ball”? may be large, small, white, black, made 
of leather or of steel, but it always has a characteristic shape. 

Some attributes of an object, denoted by a given name, occur 
always in 100% of the cases. Some others occur less frequently. Con- 
sidering each of those attributes as a component stimulus of the pat- 
tern, we may say that those component stimuli occur with different 
frequencies. We may arrange them in order of the decreasing fre- 
quency of occurrence and represent them as equi-distant points on a 
straight line, which we choose as abscissa. If we plot as ordinate the 
corresponding probability of occurrence, we shall obtain a line repre- 
senting a monotonically decreasing function. The variable along the 
abscissa may be measured in terms of the arbitrary chosen distance 
between the points. Denote the abscissa of each point by x;. If a is 
the distance between any two points, then 


a (33) 


Since each x; represents a component stimulus of the pattern, we may 
speak of “the stimulus 7;’’. 

The conditioning of the stimulus pattern S( %,, %,°+- %,) toa 
given verbal reaction R takes place gradually, as the pattern S(z,, 
%_,+++ £,) is experienced a sufficient number of times. Each of the 
component stimuli x; of the pattern S is thus conditioned to R. It 
is natural to assume that the intensity R(x;) of the reaction R to the 
stimulus x; will be the weaker, the smaller the frequency of occur- 
rence of x; in the total stimulus pattern S(2,, 22, +++ &n). 

Let us now consider a mechanism, discussed previously (Rashev- 
sky, 1936, 1938; Householder and Amelotti, 1937). Let all efferent 
conditioned fibers send inhibitory fibers to the afferent centers of each 
stimulus «;. Then a weaker stimulus x, may be completely inhibited 
by the strong reactions R(2,). Let N(#;) represent the frequency of 
occurrence of «;. Then the average intensity J(x;) of the stimulus 
x; is proportional to N(x;) and equal to b,N(x;). The intensity R (2;) 
of the reaction R in the absence of inhibiting fibers is proportional to 
I(«;) and is equal to 
R(x) = b1 (x). (34) 


But since each x; center receives inhibitory fibers from all conditioned 
pathways x; — R, therefore, the actual excitation I(x;) is given by 


I(%:) =b.N (ai) —AZ R(x), (35) 


A being a constant. 
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When the number of component stimuli is very large, we can re- 
place the sum by an integral, and thus obtain 


I(x) =b,N (e) — | R(x)dz. (36) 
0 
Introducing expression (34) into (36) we obtain 
I(x) =b,N (x) —Ab f I(x) dx. (37) 
0 


The upper limit is determined by considerations that those x’s, which 
are completely inhibited, do not produce any R(x) and do not there- 
fore contribute anything to the inhibition of others (Rashevsky, 1938, 
chap. xxviii). 

As an illustration we shall discuss here the case when 


Nir): =e, (38) 


- A component stimulus corresponding to x = 0 occurs then every 
time, the probability of its occurrence being 1. This assumption sim- 
plifies the following equations considerably, but does not introduce 
any limitation in principle. We then have (Rashevsky, 1938, pp. 282- 
283) 


I(x) = b, (e6 — er), (39) 
where « is the root of the equation oa 
oe (40) 
ba 
It is readily seen (Rashevsky, 1938, p. 283) that « decreases with in- 


creasing /. 

We can find two approximate solutions of equation (40), one for 
the case that f is very large, 6/bA > 0; and another for very small 
values of . 

In the first case we may neglect in equation (40) 6x and 1 com- 
pared with e-®*. We then have 


efk = —., (41) 
or 
E ] >0 (42) 
=— log — : 
K B ea 


In the second case we expand e* , omitting terms of third and higher 
powers. This gives 
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= ee (43) 
bpA 


In any case, for very large values of 1, that is, for very strong 
inhibition, « will be very small, and only those component stimuli will 
remain uninhibited which have a very small x, that’ is, which occurs 
either always or almost always. 

With the above mechanism only such components of S will be 
associated with the reaction R , which occur either always or almost 
always. Such components may be termed “essential’’, the remaining 
ones, “accidental”. We have here a rudimentary mechanism of ab- 
straction. The verbal reaction R is now associated not with the whole 
accidental complex of component stimuli, entering into S , but only 
with the essential ones. Thus a mathematical biophysicist considers 
as the only essential feature of the cell its metabolism, this being the 
only feature which invariably occurs in all cells (Rashevsky, 1938, 
1940). An individual with a large 4 and hence a small x will abstract 
only such features which occur always. Individuals with smaller 4 
may, at least off hand, include amongst “essential” features of an ob- 
ject, which occur very frequently, but not always. Individuals with 
very small 4 and large « , will show inability to abstract, and will con- 
sider as essential components even such which occur rather unfre- 
quently. 

If 2 is constant, then « is constant, too. Any component whose 
frequency of occurrence y(<1) is greater than e-®* would be consid- 
ered as essential, all others, as accidental. But 4, which is a meas- 
ure of the inhibitory effect, may fluctuate itself due to statistical 
physicochemical fluctuations in the organism. Then « will fluctuate 
also. Under those conditions a component stimulus with a frequency 
y of occurrence may sometimes be considered as essential, sometimes 
as not. We may now compute the probability P. for a stimulus of a 
given y to be considered essential, or the probability P, = 1— P. of 
it being considered unessential. 

Let the probability of 2 having a value between 4 and 4 + di be 
p(i)da, with p(0) = p(«o) = 0, and consider the most likely case 
that p(4) has one and only one maximum at 4 = 4,. The probability 
P:(x)d« of x having a value between «x and x + dx is obtained by sub- 
stituting into p(A) dd the expression of 4 and dd in terms of « and dx, 
obtained from either equation (42) or (43). We shall use here as an 
illustration equation (43). It gives 


2... 5h See ‘a 
BB x’ bp x wa 


A= 
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Hence 


Pils) de=— »( ma — (45) 


The abscissa x of a component of frequency y is 
1 
ay ial els (46) 


The component is considered as essential everytime when x > 2,. 
Hence the probability P, is given by 


oo 


dawlos Date \d: 
P.= Di (x) dk = — — = 
i: b B eae A (47) 


-(1/B) logy 


Expression (47) gives us P, in terms of £8, which characterizes 
the external stimulus S, and in terms of the parameters of (A), 
which contain the individual constants, in particular, 2. If p(A) is 
some kind of an exponential, P. could not be evaluated in a closed 
form. If, however, the fluctuations of 4 are approximately symmetric 
with respect to A, , then we may put 4 = A, + 2’ and consider that the 
probability of a given i’ is represented by equation (12), in which we 
ol — = i’. We now have, by substituting in equation (48), 2’ + do 
or A: 


2 4d 
Aa di’ = — zp 
b 6 x? b BK 
Due to the approximation introduced in equation (12), we now have 
eS 1 for ty < do — Va/2 and P. = 0 for 2, > & +-Va/2. For 
dg = VWa/2 < Xy < Ay + Va/2 we have 


se (48) 


EO cs ae 2 1 
b B Eh fe ey: 


—(3/B) logy 


Pe= (49) 


With p(2’) given by equation (12), P. can be expressed in closed form. 

Equation (47) or (49) can be compared with experiments. The 
quantity 6 can be controlled in principle. Thus we may show to a sub- 
ject a large number of times the same object with different accidental 
characteristics. In a selected number of cases, for instance, the ob- 
ject will be green; in another selected number of cases it will be red; 
in some cases it will be made of wood; in some of metal, etc. We thus 
can artificially “construct” the function e**, or for that matter, any 
function N(x). Expressions for P, and P. may be derived in prin- 
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ciple for any N(x). After having presented the object a large num- 
ber of times, we may ask the individual whether a given property, 
say, green, is an essential characteristic of the object or not. By re- 
cording the frequency of different replies, we may compare equations 
(47) or (49) or any other corresponding equation with experiment. 
We thus have a method of measuring quantitatively the ability to ab- 
stract. From such measurements some individual parameters of the 
subject may be calculated. 

In the following we shall consider the case of such a large 4 that 
only those-component stimuli which occur always in S remain ex- 
cited. 

Let us denote by S’ a stimulus pattern which consists of any 
stimuli lying in the interval (0,«). Denote by S” any stimulus pat- 
tern composed of stimuli lying outside of the interval (0,«). We see 
that while S’ is always contained in S, S” is only sometimes contained 
in S. In other words, we may say that “all S’s imply S’”, but only 
“some S’s imply S” ”’. 

Suppose that each peripheral afferent neuron sends off two 
branches, one to a goup A of neurons, another to a group S. Let A 
and S be identical in structure except for one thing. While each neu- 
ron of A receives from each corresponding efferent an inhibitory fi- 
ber, such as discussed above, the neurons of S do not receive such in- 
hibitory fibers. Let, however, each neuron of group A send to its 
corresponding neuron of group S a strongly inhibiting fiber. Accord- 
ing to what we have said above, when a stimulus pattern S is pre- 
sented, the individual neurons in A which correspond to any stimulus 
S” will be completely inhibited, while those corresponding to stimuli 
S’ will be excited. In S, therefore, all neurons corresponding to S’ 
will be inhibited, while those corresponding to S” will be excited. Let 
nerve pathways lead from both A and S to different conditioning cen- 
ters; and let each pair of pathways leading from A and S to any 
given conditioning center cross-inhibit each other so that only one of 
them can be conducting, the one which is more strongly excited. Let 
the neurobiophysical process corresponding to our experiencing of 
the concept “all S’s” consist in a stronger excitation of the pathway 
leading from A to any conditioning center; while the process corre- 
sponding to our experiencing of the concept “some S’s’’ consists in a 
stronger excitation of the pathways leading from S. 

If now any of the stimulus patterns S’ or any of their parts are 
associated through the proper conditioning center to any other stimu- 
lus pattern S,, then the neurobiophysical process corresponding to 
the concept “all S’s” will evoke also the stimulus S,. If any of the 
patterns S” or any of their parts are associated with any other stimu- 
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lus S. , then the neurobiophysical process corresponding to “‘some S’s” 
will evoke S.. We thus have a neurobiophysical mechanism corre- 
sponding to the notions, “all S’s imply S,” and “some S’s imply S.”. 

If a stimulus pattern S as such is presented, without specifica- 
tions such as “all” or “some”, then, according to the above, the A and 
S centers will be equally excited, and due tc mutual inhibition, none 
of the pathways from either A or S will transmit. This would mean 
simply that “S” cannot imply anything. This is not as paradoxical 
as it sounds. When we say “S implies S,” we really mean “all S’s 
imply S,”, the conventional rules of the language permitting the omis- 
sion of “all”. In other words, due to proper conditioning, omission of 
any specification results in an increased stimulation of A. 

The stimuius S,, produced indirectly, or “by implication”, itself 

will excite corresponding groups A, and S,; through any of those a 
stimulus S; may be “implied”, etc. 


ALL S;s ALL S2's ALL ao 


eee Nee 
\ /\ 
ees SS SS 28 
SOME S,'5 SOME S255 SOME “ys 
FIGURE 5 


Now consider such a chain of double pathways in which every 
preceding stimuli implies every following, as shown in Figure 5. The 
A and S pathways are symbolized by single fibers. Let each A fiber 
excite strongly every other A fiber, with which it synapses. But let 
an S fiber be only weakly excited, or not excited at all, by either an 
A or an S fiber. This can be achieved by making all S fibers have 
rather high thresholds, or by having two synapses excite every A 
fiber, but only one, every S fiber (Figure 5). In that case we have the 
following scheme: 


All S,’simply S.; all S.’simply S;; all S,’s imply S;. 


Some eon “ac Peek : “ec “ec “ec “ce some “cc 6c “ec 
All ie See as ORIG e - « «~~ no implication 
Some “é “cc “é “cc “cc “ce “é 6c“ i3 


The above scheme represents correctly some simple modes of 
logical reasoning. 

We may consider a chain consisting of any number n of such 
“elements” or “links” as represented on Figure 5. A functioning of 
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such a chain will correspond to a reasoning by a simple successive 
chain of syllogisms. 

In the further development of the theory it will be necessary to 
suggest a mechanism for “No S implies S,” and for “Some S’s do not 
imply S,”. Instead of two parallel chains of cross-inhibitory fibers, 
we shall have four chains. The treatment may have to follow the 
way indicated for the theory of discrimination between several stim- 
uli by H. D. Landahl (1938). The problem may also offer a test for 
the heuristic value of the method of W. S. McCulloch and W. Pitts 
(1948). Those authors show that by applying logical calculus, it is 
possible to construct any complicated network having given prop- 
erties. One could attempt to construct by the method of McCulloch 
and Pitts a network that would represent all modes of logical reason- 
ing, and then apply the usual methods of mathematical biophysics to 
derive some quantitative relations between different manifestations 
of the processes of logical thinking (Landahl, McCulloch, and Pitts, 
1943). 

Leaving all such developments for the future, we shall derive 
here some quantitative relations for the simple and incomplete scheme 
suggested above. Those relatiors are perhaps to be considered now 
as illustrations only. 

Each of the links of the chain, represented in Figure 5, is essen- 
tially a circuit similar to the one represented on Figure 1 and has 
accordingly a characteristic threshold h. We shall now consider the 
case when all h’s are the same. The differences in the intensities of 
excitation of the A’s and S’s correspond to the difference 4 = e, — e, 
in Figure 1. We shall assume that e, — «. is a constant quantity, the 
same for all links. This will happen if the afferent ends of the fibers 
in each link will reach their saturation values #, for rather small 
values of the stimulus (Rashevsky, 1938, chap. xxii). The synapses 
at which «, and s are considered shall be referred to as synapses 
8, and s. correspondingly. The intensities of excitation of the fibers 
going from those synapses will be denoted by E, and E.. 

Let each link send inhibitory fibers to all others, those inhibitory 
fibers rising the thresholds h. If n is the number of active links in 
the chain, the threshold of each will be given by 


h=h,+a(n—1), (50) 


where h, and a are constants. 


We may also consider the following interactions of the different 
links. 


Let the relation between « — h and E at the synapses s, and s, 
be not linear, but represented by a curve which is convex upward 
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(Rashevsky, 1938, chap. xxii). Then an equal increase in «, and é 
will have the same effect upon E, — E. as if ¢, — e, were decreased. 
If each link is connected with each other by excitatory fibers acting 
on the synapses s, and s., then the larger n, the larger «, and «, : 
and the smaller the difference FE, — F,. But it is that latter quantity 
rather than e, — e that determines 4. Only when ¢ « E-can A be 
substituted for H, — E.. Hence now the larger n, the smaller the 
“effective” 4. Within a limited range things happen formally as if 
A decreases with n , and we may put approximately 


Ae As 0th — 1). (51) 


We have the following expressions for the probabilities p., pa; 
and », of a correct, doubtful, or wrong response for each separate 
link of the chain (Landahl, 1938): 


ForA>h Fora <i 
Bo i se : i Se hira) : (52) 
4 —e™ sinh kh; a= 1—e*™ cosh kA; (53) 
Dw = 3e*O); Dw = gore, (54) 


We now may compute the probabilities P., Pa, and Py of a cor- 
rect, doubtful, or wrong response for the whole chain of n elements. 

Barring for the present the relatively unlikely case that several 
wrong conclusions combined give a correct one, we find that a correct 
total response requires a simultaneous correct response in all links. 


Hence 


A wrong response in the whole chain is obtained when either one — 
link gives a wrong response, all the others giving a correct one, or 
any two links give a wrong response, all others a correct one, etc. 
But if a single link in a chain does not respond at all, the whole chain 
will give no response. Hence, by using an argument similar to that 
which led us to equation (31) and remembering that the constants 
of all links are now assumed to be equal, we find 

n(n—1) 
P a NPP ™* ae eo a De De 
n(n —1) (n—2) 
imine tate e p 


Similarly, an expression for P, may be derived. We always have 
i C + P. ww +P d— Le. 


(56) 
| Pa in oe 
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Introducing expression (50) and (51) into expressions (52), 
(53), and (54) and the latter into expressions (55) and (56), we thus 
obtain expressions for P., P,,, and P, in terms of the number n of in- 
ferences in the whole chain of reasoning, as well as of the “individual” 
constants 4, h, a,b, and k&. Such expressions may be verified ex- 
perimentally, and from the comparison of the theory with experiment, 
some of the individual constants may be computed. 

Similar expressions may be obtained for a more general case 
when the process of reasoning is not necessarily constituted of a sim- 
ple chain of consecutive inferences. Expressions (50), (51), and (52) 
—(54) still hold for each individual inference. In general expression 
(55) for P. also still holds, but the probabilities P,- and Pz must be 
now computed in a different way depending on the “structure” of the 
logical problem. Now, for instance, some inferences of the pattern 
may give no response, yet the pattern as a whole still results in a 
total inference. For each type of logical problem we may, however, 
compute the values of P., P,, and P;, and again we may compare 
our theoretical results with experimental data. 

It must be emphasized that the derivations of equation (55) and 
(56), or any other similar expressions, is entirely independent of the 
assumed mechanism of the individual link, which provides for a sim- 
ple inference. The same method for deriving the probability of cor- 
rect, wrong, or no response in a chain of syllogism may be applied to 
any other mechanism, which is assumed for the elementary act. 

This work was aided in part by a grant from the Dr. Wailace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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In a series of recently published papers E. B. Wilson and his col- 
laborators have investigated some of the mathematical questions connect- 
ed with epidemics. A differential difference equation was set up and 
various properties of its solution were obtained. This equation implies 
that the epidemic has an equilibrium state, and it is the purpose ‘of this 
note to show that this equilibrium is stable in the sense that any small 
deviation from it will tend to zero. 


If we let S(t) denote the number of susceptible persons present 
at time t, A denote the rate of recruitment of new susceptibles, +r 
denote the time after infection when an infected person becomes in- 
fectious, « denote the duration of the disease, and 7 denote the num- 
ber of susceptibles being infected per infectious person per suscep- 
tible person per unit of time, then it has been shown (Wilson and 
Worcester, 1944) under the assumption that A ,7,o, and 7 are con- 
stant that 

A —dS/dt=r§(t) [A o—S(t—7) + S(t—7r—«)]. (1) 
The equilibrium value of S is clearly S = 1/re if Aro > 0. We shall 
explicitly make this assumption. The method of determining whether . 
or not this equilibrium is stable is standard (see, for example, Gour- 
sat, 1927). We write the equation of variation of the equation (1), 
namely 


da /dt=— Arex + o [a(t—1) —2(t—7—<«) ], (2) 
and seek to find solutions of this equation of the form x(t) = e%. 


Such a function will be a solution of equation (2), if, and only if, a 
is a solution of the secular equation 


e-oT — e-u(t+0) — (a + Are)o : (3) 


The equilibrium will be stable in case every solution a of equation (3) 
has a negative real part. This fact we shall now prove. 
Suppose first that a is a nonnegative real number. Then it fol- 


lows from the inequalities, 
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e-7 (1 —e-"7) S1— Ee" Sade < (a+ Aro)o, 


that a can not satisfy equation (3). Similarly, if a is a negative real 
number, then we have 


— eT (e- _J) <— (er — Ll) <tex= (a + Aro)c. 


It follows that equation (3) has no real roots at all. To find the com- 
plex roots, let a = b + ic, where b and ¢ are real. Then equation (3) 
is equivalent to 


eT cos ¢r — €?'T*") Cog C(7 +o) = (b+ A7ro)a, (4) 
e (749) sin e(r +o) —e sin Cr —Co. (5) 


Since there are no real roots, we have that c # 0, and we may assume 
without loss of generality that c > 0. Equation (5) may be written 
in the following form: 


f Cid — i : e(c cos cx — b sin cx) dx—=co. (6) 


We shall show that g(x) = ¢ when b 2 0, with the strict inequality 
holding at all except a finite number of points on the interval 
(7,7 +c). This will show that equation (6), and consequently equa- 
tion (5), can not be true when b2 0. 

Now g' (x) = e"[ (b? — c?) sin cx — 2be cos cx] vanishes if, and 
only if, sin. cz = +. 2b¢/ (6? +c), cos ce = = (F — ©) (6 ee 
and then g(x) = + ce™’ =c. It follows that the relative maxima of 
the function g(a) are ali less than or equal to ec. Therefore, g(x) is 
itself less than ¢ except for at most a finite number of points on the 
interval (7,7 +0). 

We have now shown that any solution (b,c) of equations (4) 
and (5) is such that b < 0. Since b is the real part of a, the equi- 
librium solution of equation (1) is stable. Moreover, since c #0, a 
small perturbation from the equilibrium will have roughly the form 
of an exponentially damped oscillation. 
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A REINTERPRETATION OF THE MATHEMATICAL 
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The fundamental equations for the interaction between neurons 
used in mathematical biophysics seem at first incompatible with the ac- 
tual neurophysiological findings on the synaptic transmission. It is 
shown, however, that those equations may be readily interpreted in 
terms of accepted neurophysiological views. What has been termed 
synapse” in mathematical biophysics must be regarded as a compli- 
cated network of internuncial neurons. it is shown that, under rather 
general conditions, the number of those interneurons will statistically 
vary with time according to the differential equation postulated for the 
excitatory and inhibitory factors. The latter are thus interpreted as the 
number of excitatory and inhibitory interneurons. 


When the fundamental postulates of the mathematical biophysics 
of the central nervous system were first developed in abstracto (Ra- 
shevsky, 1938), they were suggested, on the one hand, by the earlier 
observations of E. D. Adrian (1932) regarding the relation between 
frequency of discharge in a peripheral fiber and the intensity of the 
stimulus; on the other hand, by a possible generalization of the two- 
factor theory of peripheral excitation. The earlier abstract develop- 
ments of the theory were not too much concerned with the direct re- 
lation of the fundamental postulates to actual observations: the main 
interest centered around the proof that very complicated phenomena 
could be systematically described by a rather simple system of postu- 
lates. When, however, some of the consequences of the theory, elab- 
orated principally by A. S. Householder and H. D. Landahl, turned 
out to be in good agreement with actual observations (Rashevsky, 
1940), it became evident that the postulated equations are something 
more than mere mathematical assumptions. Yet those postulates 
seemingly did not agree with the accumulated experimental evidence 
on the interaction of two neurons at the synapse, as well as with other 
observations. It was pointed out, therefore (Rashevsky, 1940, p. 125 
ff.), that the postulated two factors, the excitatory « and the inhibi- 
tory 7, may actually be not simple physicochemical entities, but rath- 
er complex ones, and that the fundamental equations governing their 
variation with respect to time and their action may be considered as 
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convenient and useful formal postulates which eventually would have 
to be interpreted through some perhaps rather complicated mechanism. 
In the meantime, the development of the theory could go ahead, re- 
gardless of the actual true mechanism, just as the development of a 
number of formal physical theories proceeded before actual atomic 
interpretations were formed. 

Eventually the usefulness of the postulates became firmly estab- 
lished (Householder and Landahl, 1945; Rashevsky and Brown, 
1944a, b). At the same time W. S. McCulloch and W. Pitts (1943), 
using Boolean algebra, developed a theory of nervous activity which 
was based more directly on observations about synaptic transmission. 
All the successful applications of the older postulates of mathematical 
biophysics of the central nervous system remained outside of the 
scope of this new theory, although some suggestions as to how the 
apparent gap between the two approaches may be bridged were made 
by H. D. Landahl, W. S. McCulloch, and W. Pitts (1943; c.f. also 
Householder and Landahl, 1945). It is the purpose of this paper to 
show how the two approaches can be connected, and how the funda- 
mental postulates used hitherto can be reconciled with neurophysio- 
logical findings. 

The relation between frequency of discharge and intensity of a 
continuous stimulus postulated in mathematical biophysics holds only 
in a limited sense, and in a limited number of cases. For slowly 
adapting fibers, such as those of the muscle end organs (Fulton, 1943; 
Mathews, 1931a, b; 1933), the postulated relations hold for a constant 
stimulus except for the short period of the initial volleys of higher 
frequency. They also may hold in some cases of visual stimuli (Hart- 
line and Graham, 1932). For the fibers of the acoustic nerve, the rela- 
tion between the intensity of the stimulus and frequency of response 
is very similar to that postulated in mathematical biophysics, but the 
thing is complicated by the fact that the frequency of the sound en- 
ters as a parameter (Galambos and Davis, 1943). On the other hand, 
for rapidly adapting fibers such as those conveying the sensation of 
touch (Cattell and Hoagland, 1931), the postulated relations do not 
hold even with remote approximation. 

Thus while for slowly adapting fibers the postulated relation be- 
tween frequency of discharge and intensity of stimulation (Rashev- 
sky, 1938), ° 


S 1 
»=ah log or ae ark ek Ba (1) 


may hold fairly well for constant S, they certainly do not hold for 
sufficiently rapidly varying values of S. Formally, however, this may 
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be remedied by considering the rate of change S = dS/dt itself as a 
stimulus and considering a fiber as reacting both to the stimulus S 
and dS/dt. Very rapidly adapting fibers like those of the tactile re- 
ceptors, which do not discharge at all for continuous stimuli (Cattell 


and Hoagland, 1931), may be considered as responding to S only. It 
can be shown that these formal assumptions are made plausible bio- 
physically by considerations of the two-factor theory of peripheral 
nerve excitation. Such considerations also throw light on the facts 
mentioned in the preceding paragraph. 

An actual stimulus usually involves a very large number of fi- 
bers, this number increasing with the intensity due to the distribu- 
tion of thresholds. Therefore, even if we have a number of perfectly 
adapted fibers, each of which responds to a continuous stimulus with 
a frequency independent of the intensity S, still the total number »7 
of impulses per unit time for the whole set of fibers will increase with 
S, the relation between »; and S being given for physiological rea- 
sons by a curve convex upward, except perhaps in the neighborhood 
of the threshold. Hence, if the actual units which correspond to the 
“fiber” of the mathematical biophysicist are not fibers but large 
groups of fibers, then for such groups the relations (1) will still hold 
with sufficient approximation, h denoting now the lowest threshold 
of the group. 

The problem of synaptic transmission is more complicated. The 
synaptic delay between two actual neurons is of the order of 0.5 ms, 
and practically independent of the intensity of the stimulus (Lorente 
de N6, 1935a, b, c, 1938a, b, c). The frequency of a presynaptic 
discharge is not modified at the synapse, at least in sympathetic gang- 
lia (Bronk, 1939). On the other hand, the time constants for the ex- 
citatory and inhibitory factors, calculated by H. D. Landahl (1939; 
also Householder and Landahl, 1945) from comparison of theory with 
experiment are found to be sometimes of the order of seconds. Ap- 
parently, the abstract “synapse” of the mathematical biophysicist 
must involve a very large number of actual neurons. In the follow- 
ing, for precision, wherever we use the word synapse, we shall mean 
the actual synapse of the neurophysiologist. The word “synapse” in 
quotation marks will refer to the abstract concept of the mathematical 
biophysicist. 

Consider either a single very slowly adapting fiber, or a group of 
parallel fibers J (Figure 1) stimulated with an intensity S. The total 
average number of impulses per unit time in that fiber, or group of 
fibers, will then be a monotonically increasing function of S, de- 
scribed approximately by one of equations (1). Let each fiber of the 
group send off a large number of collaterals, each forming a chain of 
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neurons I’ , as indicated in Figure 1. Each chain may in general con- 
tain a different number of neurons. Let the axons of the last neurons 
of each of those chains end on a surface P, and let the collaterals of 
all the fibers J be thoroughly “mixed” in the surface P so that each 
element of area AP receives collaterals from all the fibers J. Due to 
slight variations in the individual synaptic delays of the different neu- 
rons, even a regular sequence of impulses in a fiber J will not result 
in a set of synchronous impulses arriving at the surface P. This will 
hold a fortiori for the impulses arriving at P from different fibers I 


FIGURE 1 
For simplicity, only two fibers of pathway I and their 
branches are shown in the drawing. For the same reason, 
only three circuits J” are shown. 


if the frequencies of discharges of the latter are statistically indepen- 
dent. Hence the surface P will receive a large number of impulses per 
unit time, this number being constant only on the average. But this 
average number »y will be proportional to the total average frequency 
vp in the group of fibers J . 

Let each axon of the system J’ form in the surface P a synapse 
with a neuron I”, which forms a closed chain of two or more neurons. 

Let each axon of I' excite the corresponding neuron I” just above 
the threshold so that relatively slight accidental fluctuations of the 
thresholds of I” may make them occasionally unexcitable by I’. Let « 
denote the probability that an impulse arriving at P excites the corre- 
sponding 1”. 

Consider the case in which all the neurons 7” are unexcited be- 
fore the application of the stimulus S. After S is applied, a circuit 
connected to J", is thrown into a permanently excited state every time 
an impulse arriving at P excites a neuron I”. If N, denotes the total 
number of neurons J” available at P , and N; denotes the number al- 
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ready excited, then the number of circuits thrown into an excited state 
during the time dt is dN; = a »(No — N;)dt. If, for reasons which 
will be explained presently, N; always remains much less than N,, in 
other words, when the total reserve of available neurons No is very 
large, then we have approximately 


dN; - 

—=arN,. (2) 
d 

But » is proportional to v7, and the latter is taken to be a measure of 
the intensity E of excitation of the whole pathway I. Hence, denot- 
ing by A a constant: 


dN; 
dt 


Each of the circuits, once excited, would remain so indefinitely 
if everything within it and its environment were constant. But such 
a constancy as a rule never occurs in biological systems. Thresholds 
fluctuate even in controlled experiments (Pecher, 1939; Landahl, 
1941). It may be that the number of terminal bulbs exciting each 
neuron in the circuit is so much larger than the necessary minimum 
one (Lorente de N6, 1939, p. 405) that every neuron receives a highly 
superthreshold excitation. A slight variation in the threshold in such 
a case will not affect the reverberation of the circuit. But if the num- 
ber of the terminal bulbs is just sufficient to excite, then a slight vari- 
ation in threshold, or an accidental failure of only a few terminal 
bulbs to function will break the circuit and bring it back into the un- 
excited state. It may be worth calling attention to the fact that at 
ordinary values of the pH there is about one hydrogen icn per termi- 
nal bulb. The number of molecules of other substances may be also 
very smal]. The physicochemical conditions of a single terminal bulb 
must therefore fluctuate very strongly. 

We thus see that depending on the number of the terminal bulbs, 
the circuits will possess a greater or lesser degree of instability. Due 
to accidental fluctuations, there will always be a probability of a num- 
ber of excited circuits to be spontaneously broken during the interval 
of time dt. The probability for a given circuit to fail within the in- 
terval dt is constant and equal to a dt, a being a coefficient of propor- 
tionality. Hence, the natural “rate of decay” of the excited circuits 
will be proportional to their total number N;. Therefore, when due 
to impulses arriving from I’, a certain number of circuits is thrown 
per unit time into the excited state, the rate of change of N; is given 


by 


= AE. (3) 
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dN; 
dt 


But this is formally identical with the differential equation for e or j. 

If the re-entrant part of each circuit sends off excitatory fibers, 
as shown in Figure 1, which converge upon a group of fibers JJ , then 
the total intensity of excitation of JJ will be at any moment a linear 
function of the number N; of excited circuits. On the other hand, if 
the circuits send off inhibitory fibers to IJ, we have an inhibition pro- 
portional-to N;. In view of D. Lloyd’s (1941) work, we assume with 
W. S. McCulloch and W. Pitts (1943) the existence of specific inhibi- 
tory fibers, rather then explain inhibition by the action of special in- 
ternuncials. This assumption is, however, irrelevant for the general 
argument. 

Equation (4) holds only approximately when N; << No, or 
which is the same, when AH /a << N,. For very high intensities E 
of excitation, equation (4) will have to be modified. The fact that 
equations of the form (4) can be used successfully for ¢ and 7 within 
a rather wide range of EF indicates that N, is very large. 

If no other limitations, except the obvious one, N; < N,, are im- 
posed upon N;, then we have instead of equation (4) the following: 


dN; 
dt 


= AE — aN. (4) 


=ayN,— (av +a)Nj. (5) 


Since y « FE, so that a y= cE , we may write the solution of equa- 
tion (5) for a constant EF and for N;(0) =0 as follows: 


= cN, E 
cH +a 


i [1 pees e- (ca) ] 5 (6) 

The time constant now depends on F and hence on the intensity 
S of the stimulus. Equation (6) reduces to the usual form for « 
(or 7) when cH << a;cN, =A. It would be of interest to look for 
possible indications of an increase of the time constants for very 
strong stimuli. This would show, for instance, on reaction times, 
which should decrease more rapidly for very high values of S than 
predicted by the usual theory (Landahl, 1938; Householder and Lan- 
dahl, 1945). 


——— 


FIGURE 2 
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: Instead of assuming that neurons I” form parts of closed circuits 

[interneurons C of R. Lorente de N6é (1938c, p. 210; 1939, p, 427) ], 
we may consider them as having a structure shown on Figure 2 [Type 
M of R. Lorente de N6 (1938c 210; 
, p. 210; 1939, p. 427) ]. Such a group, 
if excited at a time t, will send a regular train of impulses to a neu- 
ron of IJ on which it converges, for a time equal to ro , where r is the 
number of branches and o the synaptic delay. Such a group has al- 
ways a finite life span equal to 7c. If, however, the probability of the 
spontaneous failure of any one synapse of the group is large enough, 
the “natural” life span may practically never be reached. If the “av- 
erage” life span of such a group, determined by the probability of 
failure of a synapse, is much smaller than the “natural” life span, 
then again the rate of decay of such excited groups will be propor- 
tional to their number. 

We thus are led to interpret ¢« and 7 as a measure of the number 
of excited groups of interneurons of a certain type. Since the sta- 
bility of such groups may vary within a very wide range, we obtain 
a very wide range of variation for the “synaptic” delay referred to 
in mathematical biophysics. The “synapse” of the mathematical bio- 
physicist includes all interneurons between J and IJ (Figure 1). In 
some cases it may include a large portion of the brain. This throws 
light on the puzzling situation that reaction-times are correctly de- 
scribed quantitatively by the mathematical biophysics of a two-neuron 
are with a single “synapse”, although actually a large number of 
synapses is involved (Rashevsky, 1940). 

In the light of all the above, it is suggested that in mathematical 
biophysics the word “connection” be substituted for “synapse”, and 
the word “pathway” for “fiber”. The word “neuron” should be re- 
tained only in its exact anatomical sense. “Group of neurons” or 
“macroneuron” should be used otherwise. 

In all cases, when we deal with very large numbers of afferent 
and efferent neurons, we may continue to use the fundamental equa- 
tions of mathematical biophysics, with possible improvements that 
may be suggested by the present interpretation. Whenever we deal 
with phenomena involving relatively few neurons the method of W. 
S. McCulloch and W. Pitts will be used. The situation reminds us of 
the parallel use of thermodynamics and kinetic theory in physics, de- 
pending on the nature of the problem treated. 

In particular, in discussing closed circuits, we must keep in mind 
the difference between the “microscopic” and macroscopic” circuits, 
or “circuits of individual neurons” and “circuits of pathways”. A 
single neuron cannot form a self-exciting circuit because the excita- 
tion will fall within its own refractory phase. But a group of neu- 
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rons, or a “macroneuron” may well be circuited in itself (House- 
holder and Landahl, 1939). It is likely that/ such circuits as were 
designated by Eng in a previous paper (Rashevsky, 1945) are to be 
conceived as “microcircuits”, while closed circuits discussed in the 
same paper and assumed to connect different regions of the brain 
(loc. cit., Figure 7) are to be considered as circuits of pathways. 
One interesting consequence may be noted. If e and 7 measure 
the number of excited groups of interneurons, then the statistical 
fluctuations of e and j in the stationary state, such as introduced by 
H. D. Landahl (1938, also Householder and Landahl, 1945) follow im- 
mediately from the above picture. Both the building up of « or 7, as 
well as their decrease follow equation (2) only statistically. The de- 
cay, for instance, does not proceed exactly according to an exponen- 
tial. If m in the average number of groups decays during the time 
t, that is, 
Mm = aN ;t A (7) 


then (Bothe, 1926) the probability that in the interval t actually m 
groups will decay is given by 


mm — 

Dae eae (8) 

m! 

or for large values of ™ , 
(m—m)2 
iL ES 

Dees em , (9) 

V2am 


Denoting by 4, the average absolute deviation m — @ , we have for 
the relative variation, 


An/m =1/V/Mm . (10) 


Due to the fluctuation of the decay, the stationary values of « 
and 7 will also fluctuate even if the building-up rate were constant. 
The fluctuation of the latter will complicate the picture further. But 
in general, we shall find that the fluctuations will decrease with in- 
creasing e« and 7. Such relations may be in principle verified experi- 
mentally by comparing theoretical and experimental data for psycho- 
physical discriminations, as has been done by H. D. Landahl (1938; 
. also Householder and Landahl, 1945), but using a very wide range 
of intensities of stimuli. 

The author is indebted to Dr. H. D. Landahl for a critical dis- 
cussion of the paper and for valuable comments. 
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